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Abstract 

Brown's formulation of dynamical perfect fluids in Minkowski space-time is extended to ADM 
tetrad gravity in globally hyperbolic, asymptotically Minkowskian space-times. For the dust we 
get the Hamiltonian description in closed form in the York canonical basis, where we can separate 
the inertial gauge variables of the gravitational field in the non-Euclidean 3-spaces of global non- 
inertial frames from the physical tidal ones. After writing the Hamilton equations of the dust, we 
identify the sector of irrotational motions and the gauge fixings forcing the dust 3-spaces to coincide 
with the 3-spaces of the non-inertial frame. The role of the inertial gauge variable York time (the 
remnant of the clock synchronization gauge freedom) is emphasized. Finally the Hamiltonian Post- 
Minkowskian linearization is studied. The future application of this formalism will be the study of 
cosmological back-reaction (as an alternative to dark energy) in the York canonical basis. 
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I. INTRODUCTION 



By using Dirac's theory of constraints in a series of papers [1-4] we developed the Hamil- 
tonian formulation of ADM tetrad gravity in suitable globally hyperbolic, asymptotically 
Minkowskian space-times admitting global 3+1 splittings (i.e. global non-inertial frames) 
with dynamical matter represented by positive-energy charged scalar point particles and by 
the electro-magnetic field in the radiation gauge. By means of a Shanmugadhasan canonical 
transformation, implementing the York map [5] , we found a York canonical basis adapted to 
ten of the fourteen first class constraints of tetrad gravity and we could identify the gauge 
variables (the inertial effects) and the physical degrees of freedom (the tidal effects) of the 
gravitational field in global non-inertial frames (the only ones allowed by the equivalence 
principle). 

In Refs.[3, 4] we also developed the Hamiltonian Post-Minkowskian (HPM) linearization 
in a family of (non-harmonic) 3-orthogonal Schwinger time gauges. In particular it turns 
out [4] that at least part of the astrophysical and cosmological dark matter could arise as 
a relativistic inertial effect determined by the trace of the extrinsic curvature (York time) 
of the instantaneous 3-spaces as 3-manifolds embedded into space-time. This inertial gauge 
variable is the general relativistic remnant of the special relativistic gauge freedom in clock 
synchronization as a definition of instantaneous 3-space. In canonical general relativity the 
3-spaces are dynamically determined modulo this gauge freedom [6]. 

With the exception of the three gauge variables describing the freedom in the choice of the 
3-coordinates inside the instantaneous 3-spaces of the non-inertial frame all the other inertial 
and tidal effects are 3-scalar fields on the 3-spaces in this Hamiltonian formulation. The 
3-scalar nature of the tidal variables and of many of the inertial ones opens the possibility 
to study the averaging over 3- volumes of most of the Hamilton equations (and not of only a 
small subset as in the standard approaches), extending the results of Buchert approach [7] 
to cosmology in which the back-reaction arising from the volume-average of the non-linear 
Einstein equations is proposed as an alternative to dark energy. It will be interesting to see 
which is the role of the York time in a framework using neither 3-spaces comoving with an 
irrotational fluid nor comoving coordinates (see Ref. [8] for the reformulation of the approach 
of Ref. [7] in arbitrary coordinate systems) . 

But before doing this attempt we need to include dynamical perfect fluids in our frame- 
work, because this is the standard matter (very often only test matter) in cosmology. More- 
over we want to avoid the use of irrotational test fluids in cosmological comoving coordinates, 
because this simplification does not allow to make explicit the gauge freedom in the York 
time. 

In this paper we define the general relativistic extension of Brown's approach [9] to 
dynamical perfect fluids in Minkowski space-time. Our starting point will be Refs.flO, 11], 
where we reformulated this approach in the framework of the inertial rest frames, a subcase 
of the general treatment of special relativistic non-inertial frames developed in Ref. [12]. 

We consider in detail dust, because, as shown in Refs.[10], only in this case (and in few 
others including the photon gas) we can get a closed form for the Hamiltonian mass density 
of the fluid 1 . 



In general the transition from velocities to momenta in this approach leads to a trascendental equation, 
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In Section II we give a review of the 3+1 point of view which allows to define global 
non-inertial frames with well defined instantaneous 3-spaces necessary for the formulation 
of the Cauchy problem of field equations. After the special relativistic formulation based on 
Ref.[12], we give its extension to the class of space-times used in ADM tetrad gravity [1, 2] 
together with the definitions of the tetrad and cotetrad fields. In both cases there are two 
congruences of time-like observers associated with each 3+1 splitting of space-time: i) the 
congruence of the Eulerian observers, whose 4-velocity is the unit normal to the 3-spaces; 
ii) the skew, in general non surface-forming, congruence, whose 4-velocity is the evolution 
vector field associated with the 3+1 splitting of space-time. 

After having clarified this kinematical scenario we review in Section III Brown's approach 
[9] to perfect fluids in Minkowski space-time in the formulation of Ref.[10] employing the 
rest-frame instant form of the dynamics of isolated systems [12]. We show that we can 
get the explicit Hamiltonian expression of the energy-momentum tensor and of the unit 
4-velocity of the fluid only for special equations of state, in particular for dust. In the rest of 
the paper only the dust will be considered. The dust is not irrotational so that its 4-velocity 
is not surface-forming. 

Then in Section IV we study the action of ADM tetrad gravity coupled to dynamical 
matter used in Refs.fl, 2] with the matter action of a dust. After having obtained the 
Hamiltonian formulation we make the canonical transformation to the York canonical basis 
with the resulting expressions for the 4-velocity, the energy-momentum tensor and the dust 
Hamilton equations. Also the expressions in the family of (non-harmonic) 3-orthogonal 
Schwinger time gauges, where the spatial metric is diagonal, is given. Finally we present 
the Hamiltonian Post-Minkowskian linearization of Refs.[3, 4] in the case of dust. 

Section V contains a study of the skew congruence of time-like observers defined by the 
unit 4-velocity of the dust (their world-lines are the flux lines). It is shown that a gauge 
fixing of the inertial shift functions of the gravitational field (it is a restriction on the 3- 
coordinates on the 3-spaces) is needed to make it to coincide with the skew congruence of 
the 3+1 splitting of space-time. After analyzing the properties of the dust flux lines, we 
introduce the Eulerian point of view for the description of the dust following Ref.[ll]. Then 
we introduce the acceleration, the expansion, the shear and the vorticity of the dust. We 
find which are the two first class Hamiltonian constraints to be added by hand to restrict 
the dust motions to the irrotational ones. These motions are described only by one pair 
of canonical variables and their 4-velocity is surface-forming: it can be made to coincide 
with the unit normal to the 3-spaces of the 3+1 splitting by adding a gauge fixing on the 
inertial lapse function (it is a restriction on the York time, i.e. on the clock synchronization 
freedom). In this gauge the irrotational dust is comoving with the 3-space and the flux 
lines are the world-lines of the Eulerian observers. We studied all these properties in detail, 
because usually these topics are not discussed in the literature. 

In the Conclusions we make some comments on the future use of these result in cosmology. 

In Appendix A there is a review of the York canonical basis of ADM tetrad gravity 
[1, 2], followed by the definition of the acceleration, expansion and shear of the irrotational 
Eulerian observers whose 4-velocity is the unit normal to the instantaneous 3-spaces. 



whose explicit solution is known in few cases. 
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In Appendix B there the discussion on the relation between Lagrangian velocity and 
Hamiltonian momenta as a function of the equation of state of the perfect fluid. 
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II. THE 3+1 POINT OF VIEW AND GLOBAL NON-INERTIAL FRAMES IN 
SPECIAL AND GENERAL RELATIVITY 



In this Section we review some aspects of the theory of global non-inertial frames, the 
only ones admitted in general relativity due to the equivalence principle. Their formulation 
in Minkowski space-time was given in Ref.[12]. Then it was extended to globally hyperbolic, 
asymptotically Minkowskian Einstein space-times without super-translations in Refs.fl, 2]. 

A. Minkowski Space-Time 

As shown in Ref.[12] we now have a metrology-oriented description of non-inertial frames 
in special relativity. This can be done with the 3+1 point of view and the use of observer- 
dependent Lorentz-scalar radar 4-coordinates. Let us give the world-line a+(r) of an arbi- 
trary time-like observer carrying a standard atomic clock: r is an arbitrary monotonically 
increasing function of the proper time of this clock. Then we give an admissible 3+1 split- 
ting of Minkowski space-time, namely a nice foliation with space-like instantaneous 3-spaces 
S r : it is the mathematical idealization of a protocol for clock synchronization (all the clocks 
in the points of £ r sign the same time of the atomic clock of the observer). On each 3- 
space S r we choose curvilinear 3-coordinates a r having the observer as origin. These are 
the radar 4-coordinates a A = (t; a r ). If x^ h-> <j a [x) is the coordinate transformation from 
the Cartesian 4-coordinates x^ of a reference inertial observer to radar coordinates, its in- 
verse o A i — y a+ = z fJ, ( y r,cr r ) defines the embedding functions z^(t, o~ r ) describing the 3-spaces 
S r as embedded 3-manifold into Minkowski space-time. From now on we shall denote the 
curvilinear 3-coordinates a r with the notation & for the sake of simplicity. 

The induced 4- metric on S r is the following functional of the embedding: VasC*", = 
[za Vnv z b]( t i +)> where = dz^/da A and 4 r/^ is the flat metric. The 4- metric 4 gAB has 

signature e (H ) with e = ± (the particle physics, e = +, and general relativity, e = — , 

conventions); the flat Minkowski metric is rj^ = e (H ). 

While the 4-vectors z^(r,a) are tangent to S T , so that the unit normal / M (r, <?) is pro- 
portional to e M a/ 3 7 [z" z% zj](r, a), we have z%(t, (?) = [N Z M + n r z^](r, a) for the so-called 
evolution 4- vector, where N(r, a) — 1 + n(r, a) = e l^ij, a) and n r (r, a) = —eg rr (T, a) = 
[ 3 g r s n s ] (r, a) are the lapse and shift functions. We also have \det 4 g\ = (1 +71)^/7; 
y^y = a/ det 3 g with 3 g rs = — e 4 g rs of positive signature. 

The conditions for having an admissible 3+1 splitting of space-time are: 

a) 1 +n(r, a) > everywhere (the instantaneous 3-spaces never intersect each other, so that 
there are no coordinate singularities like it happen with Fermi coordinates); 

b) the M0ller conditions [12], which imply 

i) e 4 g TT > 0, i.e. (1 + n) 2 > J2 r n r nT (the rotational velocity never exceeds the velocity of 
light c, so that there no coordinate singularities like it happens with the rotating disk); 

ii) e 4 g rr = — 3 g rr < (satisfied by the signature of 3 g r s), A grr 4 fi f ss _ ( 4 fi'rs) 2 > and dete A g rs = 
—det 3 g rs < (satisfied by the signature of 3 g rs ) so that det 4 gAB < (these conditions imply 
that 3 g rs has three definite positive eigenvalues A r = in the non-degenerate case without 
Killing symmetries, the only one we consider). 

In this 3+1 point of view the embedding functions z^(r : a) describe the inertial effects 
present in the given non-inertial frame. The 4-metric A gAB(j,&) is the potential for the 
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induced inertial effects. For instance the extrinsic curvature 3 K rs (r, a) = yjJT+n) ( Ur \ s + 

n s\r — d T 3 g rs )^(T,(r) of the non-Euclidean 3-spaces S r (|r denotes the covariant derivative 

in it) is one of these induced inertial effects. All these inertial effects derive from the gauge 
freedoms of clock synchronization and choice of the 3-coordinates. 

In Ref. [12] there is a complete description of isolated systems in non-inertial frames 
by means of parametrized Minkowski theories where there is a well defined action principle 
containing the embeddings z^ij, a) as Lagrangian variables and implying their gauge nature. 

In special relativity we can restrict ourselves to inertial frames and define the inertial 
rest- frame instant form of dynamics for isolated systems by choosing the 3+1 splitting cor- 
responding to the intrinsic inertial rest frame of the isolated system centered on an iner- 
tial observer: the instantaneous 3-spaces, named Wigner 3-space due to the fact that the 

3- vectors inside them are Wigner spin-1 3-vectors [12], are orthogonal to the conserved 

4- momentum P M of the configuration. In Ref. [12] there is the extension to admissible non- 
inertial rest frames, where P M is orthogonal to the asymptotic space-like hyper-planes to 
which the instantaneous 3-spaces tend at spatial infinity. This non-inertial family of 3+1 
splittings is the only one admitted by the asymptotically Minkowskian space-times covered 
by the canonical gravity formulation of Refs. [1-4] . 



B. Globally Hyperbolic, Asymptotically Minkowskian Einstein Space-times 

In the globally hyperbolic, asymptotically Minkowskian space-times without super- 
translations of Ref. [1] we can use radar 4-coordinates a A = (cr T = r; a r ), A = r,r, adapted 
to an admissible 3+1 splitting of the space-time and centered on an arbitrary time-like ob- 
server like in special relativity: they define a non-inertial frame centered on the observer. 
The absence of super-translations implies [13] that the instantaneous 3-spaces S r are or- 
thogonal to the conserved ADM 4-momentum at spatial infinity, i.e. they are non-inertial 
rest frames of the 3-universe. Therefore at spatial infinity there is an asymptotic Minkowski 
background J^-metric and there are asymptotic inertial observers to be identified with the 
fixed stars of the star catalogues. 

The 3-spaces £ r are identified by the embeddings x M = z^(r, a), but now the quanti- 
ties z^(t, a) = 9 Z q^a^ are the transition coefficients from the adapted radar 4-coordinates 
a A = (r,a r ) to world 4-coordinates a+ ( A gAB = z^ 4 g fiu z v B ). As in special relativity 
the space-like 4- vectors z^(r,a) are tangent to the 3-spaces, the unit normal to them 

is l"(r,a) = (^Z a )(t,<7) = V^^d e ^ z ° Z P z fj ( r> an d the time-like evolu- 

tion 4-vector z%(t, a) has the decomposition z^(r,a) = ^(1 + n) + n r z^(r, a), with 

i+(r, a) = (za v A j (r, <?) = (z>? / \f e A g^ u z? (r, <?) the associated unit time-like 4-vector 
(skew with respect to the foliation with 3-spaces). 

In general relativity, where the 3-spaces are dynamically determined [6] except for the 
York time 3 A" (the general relativistic remnant of the special relativistic gauge freedom in 
clock synchronization; differently from special relativity 3 K(t, a) is an independent inertial 
gauge variable and not an induced inertial effect), it is convenient to use the embedding 
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,2 m (t, a) = + e^r + e^cr r , where e A are asymptotic flat tetrads (e£ is orthogonal to the 
asymptotic flat Euclidean 3-space, being proportional to the conserved ADM 4-momentum). 

In Refs. [1-4] (see also Appendix A) we introduced ADM tetrad gravity by considering 
the ADM action as a functional of cotetrads 4 E A a \r, a) = z a (t, a) 4 E^\t, a) 2 by using the 
following decomposition of the 4-metric in radar 4-coordinates 

o (a) o (P) 

9AB - &a V(a)(P) J^B - V(a)(P) &B > 

o(/3) o(°) , , °<» 

4 4 a) =L^ w ( m )*E A =L^\ o) { m yE A +L$( m )Rl m (a (d) ) 4 E A . 

(2.1) 



o (a) 

The cotetrads 4 E A are adapted to the 3+1 splitting (the time-like adapted tetrad is 
the unit normal l A to E T ) by a point-dependent standard Lorentz boost for time-like orbits 
acting on the flat indices parametrized by </?(„) (r, a) 3 . The adapted tetrads and cotetrads 
(corresponding to the so called Schwinger time gauges) have the expression 



% ) = 7^ (1; -»(«) %)) = l A , % a) = (0; 3 e[ a) ), 

4 £a = (1 + n) (1;0) = el A , 4 E A = (n (a) ; 3 e (a)r ), (2.2) 

where 3 e^ a ^ (r, a) and 3 e( a ) r (r, a) are triads and cotriads on S T and ri( a ) (r, a) = 
(n r 3 e T ( a )J (t, <?) = (n r 3 e( a ) r j (r, a) 4 are adapted shift functions vanishing at spatial infinity. 

o A o A 

The adapted tetrads 4 E^ and triads 3 e r ^ are defined modulo SO (3) rotations 4 E^ = 

R( a )(b)(u(e)) 4 -E(b), 3 e( a) = R( a )( b )(a {e )) 3 e r {b) , where a (a )(r, <r) are three point-dependent Euler 
angles. After having chosen an arbitrary point-dependent origin a^(r,a) = 0, we arrive at 
the following adapted tetrads and cotetrads [n^ a ) = J2b n (b) R(b){a)( a {e)) ] 



2 Flat indices (a), a = o,a, are raised and lowered by the flat Minkowski metric 4 ?7( Q )(^) = e (H ). We 

define 4 f?( a )(&) = — e#(a)(6) with a positive-definite Euclidean 3-metric. 

3 As shown in Ref.[l, 13], the flat indices (a) of the adapted tetrads and cotetrads and of the triads and 
cotriads on E T transform as Wigner spin-1 indices under the point-dependent SO(3) Wigner rotations 
associated with Lorentz transformations A^ a hm(z) in the tangent plane to the space-time in the given 
point of S T . Instead the index (o) of the adapted tetrads and cotetrads is a local Lorentz scalar index. 

4 Since we use the positive-definite 3-metric <5( a )(6)i we shall use only lower flat spatial indices. Therefore 

for the cotriads we use the notation 3 er^ d = 3 e( a ) r with £( a )(6) = 3e L) 3e (b)r- 
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o A qA 



o A 



4 Z7' (\. ^ 3-r \ lA 4771 fn. 3-r 

#( o) = E {o) = - ■ - ••• ■ 



1 + n 



(l;-n(a)%)) = l A , % a) = (^e r (a) ), 



o (o) o (o) o («) o 

4 £ A = 4 £ A = (l + n)(l;0) = e/ A , 4 E A = 4 E (a)A = (n (a) ; 3 e (a)r ). (2.3) 



The future-oriented unit normal to S r and the projector on E T are 



l A 



e(l + n) (l; o), 4 ^ AB Z A Z B = e, 

e (1 + n) V r = — *— ( 1; -^ r ) = ( 1 

V J y 1 + n V ' / 1 + n V 



5 ~ n (a) 3 e( a) I , 



Wi = 



5 B -el A l B , 3 hl = 3 K = 0, 3 K = n {a) 3 el a) , 3 h r s = 5 r s , 



h TT e ^ ^ n (a)' h TT C^(a) C(a)rj ^rs e C(a)r e (a)sj 

r 3-s 
(a) e (a)- 



SfcTT = 3 fc rr = Qj S fc « = _ c 3gr 3-s 



(2.4) 



The 4-metric has the following expression (the lapse and shift function are independent 
inertial gauge variables) 

V = e [(1 + n) 2 - V s n r nj = e [(1 + n) 2 - £ n 2 a) ] = e [(1 + n) 2 - £ n 2 a) ], 



4 9rr = -en r = -en( a ) 3 e (a ) r = -en^ 3 e( a ) r , 

9rs C 9rsi 



9rs ' C(a)r &{a)s &(a)r ^(a)sj 



3^rs _ 3r.rs _ 3 r 3 s _ 3^r 3^s 

e (a) e (o) - e (a) e (a) , 



i^rs = 3^r 



4 g rs 



(1+n) 2 ' 

-6 (V - 



V = -6 



3e (a) "(a) 



3 e[a) "(a) 



(l + n) ; 



(1 + n) 5 



= -^ e \a) 3e (b) ( S (a)(b) 



(1 + n) 2 (1 + n) 2 ' 

n>( a ) ri(b) \ _ 



(1 +nf 



^3^r 3-s /r 3„rs _ r,rs 

- e e (a) e (6) (0( a )(6) - ^ ; ^J , g -h , 



(1 + n) 2 



4 ^TT ^rS 4 (^ rr ^<? TS 



3^s 



(1+n) 2 ' 



^ = 7=( 3 e) 2 , 3 e = det 3 e (a)r , 



v^(l + n) = 3 e(l + n). 



(2.5) 
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The 3-metric 3 g rs has signature (+ + +), so that we may put all the flat 3- indices down. 
We have 3 g ru3 g us = S r s , d A 3 g rs = - 3 g ru 3 g sv d A 3 g uv . 



C. The Two Congruences of Time-like Observers associated with Non-Inertial 
Frames 

Each 3+1 splitting of space-time, i.e. each global non-inertial frame, has two associated 
congruences of time-like observers: 



T,<7 



i) The congruence of the Eulerian observers with the unit normal Z M (r, a) = [z^ l A ^j (i 

to the 3-spaces as unit 4- velocity. The world-lines of these observers are the integral curves of 
the unit normal and in general are not geodesies. In adapted radar 4-coordinates the contro- 

o A o 

variant (l A (r, a), A E^{r,a)) and covariant (Ia(t, <t), a E^)a(t, &)) orthonormal tetrads car- 
ried by the Eulerian observers are given in Eqs.(2.3). 

ii) The skew congruence with unit 4- velocity i> M (r, <?) = (^z^v a ^(t, a) (in general it 

is not surface- forming, i.e. it has a non- vanishing vorticity). The observers of the skew 
congruence have the world-lines (integral curves of the 4-velocity) defined by a T = const. 
for every r, because the unit 4-velocity tangent to the flux lines Xg o (j) = z^{j,a ) is 

Vg{r) = z^{t,B ) I \J~e A g TT (r, a a ) . They carry the adapted contro- variant orthonormal 

o A 

tetrads (V^(r, a) are not tangent to the 3-spaces S T like A E^ (r, a) of Eqs.(2.3)) 

v A (r,a) = (1;0) =(r,a), 

nj a) 

The covariant version of these tetrads is 



ev A (r,a) = ( /(I + >»)»-£ nfc; 7=== -J^' 

V {a)A {r,a) = (0; 3 e (a)r )(r )( 7). (2.7) 

In each point there is a Lorentz transformation connecting these tetrads to the ones of 
the Eulerian observer present in this point. 

When there is a perfect fluid with unit time-like 4-velocity C/ A (r, <?), there is also the 
congruence of the time-like flux curves: in general it is not surface-forming and it is in- 
dependent from the previous two congruences. If (u a {t, a); U^(r, a) j is an orthonormal 
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tetrad carried by a flux line, the connection of these 4-vectors to the orthonormal tetrad of 
the Eulerian observers 5 is 



U A (r,B) = r(^ + ^/3 (a) 4 E (a) )(r,a), 

a 

/ °- A \ 

W£)(r, 3) = (f (o) l A + J2 7(a)(6) %)) (r, a), (2.8) 



(t,<t) 



with i (a) (r, <?) = ( E 6 7(a)(6) £(&)) (r, <?) and £ cd (<S (c)(d) - /3 (c) (tJ) ) 7(o)(c) 7(6)(d) 
^(o)(6)- See Section V for the gauge fixings implying U A {r, (?) ~ l A {j, <?). 

In Ref. [8] (see also Refs. [14]) the results of Ref. [7] are reformulated in arbitrary coordi- 
nates with the following change of notation 

U A = F(w) (l A + m A )=T(w) (l ; -[J2n ia) 3 e[ a) + (l + n)w r ]), 

a 

m A = E A.) 4 K) = (o; = E 0w Sg w)> 

a a 

C/ A = r(iu) + m A j = e r(tw) ^1 + n + E "(a) 3g (a) ^ W r) , 

as 

mA = C (E n (a) 3 e[a) W s;Wr^, m A l A = 0, 



as 



T(w) = 1 = 1 = y/1 + 3g"U r U a , 

a/1 - A g rs w r w s ^l- 6 g rs w r w s 

w r — e — , w — q w s . 

y^TTWTvUs 



(2.9) 



When the vorticity of the fluid vanishes, so that its 4- velocity is surface forming, there is a 
3+1 splitting of space-time determined by the irrotational fluid. While in special relativity 
we can always choose a global non-inertial frame coinciding with these 3+1 splitting, in 
general relativity we have to show that there is a gauge fixing on the inertial gauge variable 
3 K(t, (?) (the York time) allowing this identification (see Section V). 

Let us remark that Eqs.(2.8) establish the bridge between our 3+1 point of view and 
the 1+3 point of view of Refs. [15], where one describes both the gravitational field and the 
matter as seen by a generic family of observers with 4- velocity U a (t, <?). Most of the results 
in cosmology (see for instance Refs. [16]) are presented in the 1+3 framework. However, in 
the 1+3 point of view vorticity is an obstruction to formulate the Cauchy problem (3-spaces 



5 Analogous expressions hold in terms of the tetrads (2.6) of the observers of the skew congruence. See 
Eqs.(5.3) for the gauge fixings implying U a (t,<t) rts v A (r,a). 
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are not existing; each observer uses as rest frame the tangent 3-space orthogonal to the 
4-velocity) and there is no natural way to identify the inertial gravitational gauge variables 
of the Hamiltonian formalism based on Dirac's constraint theory (see also Appendix A of 
paper I for the induced treatment of the non-Hamiltonian ADM equations). 
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III. PERFECT FLUIDS IN MINKOWSKI SPACE-TIME 



In this Section we will review Brown approach [9] to perfect fluids in Minkowski space- 
time. In Ref.[9] a set of space-time scalar fields a l (x), i = 1,2,3, was interpreted as La- 
grangian (or comoving) coordinates for the fluid. They label the fluid flow lines (physically 
determined by the average particle motions) passing through the points inside the bound- 
ary. In an allowed 3+1 splitting of Minkowski space-time, we use the three scalar fields 
o?(r, cr) = a l (z(T, a)) for the Lagrangian coordinates of the fluid in the instantaneous 3- 
space E T . Now either the a l (r, <r)'s have a compact boundary V a {r) C £ r or have bound- 
ary conditions at spatial infinity. For each value of r, one could invert a 1 = a l (T,a) to 
a r = a r {T,a % ) and use the a*'s as a special coordinate system on S T inside the support 
V a (r) C S T : z"(t, <?(t, a')) = z^{r, a'). 

The unit time-like 4-velocity field U a (t,cj) of the fluid is a derived function of the La- 
grangian coordinates of the fluid. In this way, as shown in Ref.[9], one has not to add new 
quantities (like the Clebsch potentials) with Lagrange multipliers to the action to implement 
the particle number conservation and the absence of entropy exchange between neighboring 
flow lines: they are automatically satisfied as a consequence of the comoving nature of these 
Lagrangian coordinates, which implies 

(u A d A of}{T,S) =0. (3.1) 



Instead of the 4-velocity it is convenient to introduce a material current [9, 10] 
J A (a l (T, a)) = y/\ 4 g(r, <?)| n(r, o) U A (r, a), (^J A 8a a*) (t, a) = 0, where n(r, a) is the parti- 
cle number density. This material current has the following form 6 



J A (a l (T,a)) = [(l + n)^nU A ](r,a), [J T d T a 1 + J r d r a l ](r, a) = 0, 

J T {a\T, a)) = [e ruv d r a 1 d u a 2 d v a 3 ](T,a) = 

= i e ruv e ijk [d r a 1 d u a j d v a k ] (r, a) = det (d r a 1 ) (r, a) , 

J r \a\r,<j)) = -±e™e ijk [d T a i d u a i d v a k \(T,ff), 

{ U ] d T a l = -^d r a l = -^-d r a\ (3.2) 

with Af = f v ( T }d 3 <j J T (a l (r,a)) giving the conserved particle number and with 
f v ( r ) d 3 a(s J t )(t, a) giving the conserved entropy per particle. 

This shows that the fluid flow lines, whose tangent vector field is the fluid 4-velocity time- 
like vector field U A , are identified by a 1 = const, and that the particle number conservation 
is automatic. Moreover, if the entropy for particle is a function only of the fluid Lagrangian 



6 Here J T has the opposite sign with respect to Ref.[10]. We have e ruv e^k d m a 1 d u a 3 d v a k = 
2 5 r m det {d u a k ) = 2 J T S r m and e ruv e ijk d r a h d u a 3 d v a k = 2 J T 5?. 
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coordinates, s = s(a l ), the assumed form of J A also implies automatically the absence of 
entropy exchange between neighboring flow lines, 8a (n sU A ) = 0. Since U A dAs(a % ) = 0, 
the perfect fluid is locally adiabatic; instead for an isentropic fluid we have 8as = 0, namely 
s = const.. 

Since U A is a unit 4-vector, we have the following expression for the particle number 
density and for the unit 4-velocity in terms of the material current 



^ \A ( ~ Ve*g AB (T,a) J A (a^r,a)) jggrgg 
n(r,a) = — : — — (t, er) = r . — ; — -r-; 



1 /IT „ , J r + n r J T J s + n s J T . . 

(J T ) 2 - 3 9rs— r- — (r, a), 



a/7(t, a) V 1 + n 1 + n 

jjAr^. &) = J A (a l (T,a)) 

T/izs unii ^-velocity is in general different from the unit normal l A {r,a) to the 3-space S r 
and therefore in general it is not surface forming. 

The action given in Ref.[9] for an (isentropic if s = const.) perfect fluid with equation 
of state p = p(n, s) 7 , can be reformulated as a parametrized Minkowski theory [10, 12] if 
written in the form 



' fluid 



drd 3 a L^gAB^r, a),a\r, a)) = 



drd 3 a((l + n) ^)(r, a) p 



\J(a*(r,a) 



((1 + n) y/T)(T,ff)' 



s(a 



-J drd 3 a((l + n) ^f)(r,a) 
1 



P 



:j t ) 2 - 3 9u 



ju + n u Jr Jv + n v Jr 



1 + n 



1 + n 



(r,<r), s(o! i (r,a))). 



(3.4) 



Here is a functional of the fluid and of the embedding z^{r,cr). See Ref.[12] for the 

general treatment of the parametrized Minkowski theories, for the gauge equivalence of the 
descriptions in different admissible 3+1 splittings of Minkowski space-time and for rest-frame 
instant form description. Here we only emphasize the properties of the fluid developed in 
Ref. [10] in view of the extension to general relativity. 

The pressure of the fluid is 



/- \ ~dp(n,s) 
p{n,s) = n — — — \ s -p(n,s). 



(3.5) 



7 As shown in Ref. [9], in the case that the entropy has the form s = s(a l ) the resulting Euler-Lagrange 
equations imply the standard Euler equations obtained from the conservation of the energy-momentum 
tensor. See also Ref. [10]. 
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The canonical momenta of the fluid are 



tt/^x 6 s fimd Y r (r,a) f dp(h,s) \ l 

lk{r,(T) = -— — = — —l ri {T,cr)[ — )\n=x{T,(j), 

dd T a % (T,a) X(t,<j) \ on J 

with the notations 

JU _|_ „U JT 

X=\J\ = y/{jry - 3^ w y«y« = ^ fi, Y u = — ±- , 

1 + n 

7_ 3 SMt) Q jr Q fc <T- Q j 3 TT 

ri 2 u ^ ' " n ^ rn ' 

Y r = u i (r- 1 y r -7 — ^ = — V s ^^n^ — ^ = y r (x,nA 

(3.6) 

The matrix Th is equal to E s 3 g rs (ad J'is) , where adJi r = J T J^ 1 is the adjoint matrix 
of the Jacobian J = (Ji r = d r a 1 ) of the transformation from the Lagrangian coordinates 
cz 1 (t, a) to the Eulerian ones a on S r (see Ref.[ll] and Section III). We use the notations 
J T = (detj), E, ^rXr 1 = 5 rs , E r ^Jr/ = 5 tv E r (T' 1 ) ir T rj = E. % t (T~ l T = 5 s r . 

The main problem of this approach is to invert these equations to get the velocities 
d T a l (T,a) in terms of the momenta Hi(r, a), i.e. to get the unit J^-velocity U A {r,a) and 
then the energy-momentum tensor and the Hamiltonian in terms of the momenta. 

We have to find the solution X of the equation X 2 + 3 g rs Y r (X,Ui) Y s (X,Ui) = (J T f 
with Y r (X, rij) given by the last line of Eq.(3.6). This equation can be rewritten in the 
following form 



X 2 



JT ) 2 + A 2 {^^y Z \ fi= ^]=(r)\ A 2 = J2 3 9 rS d r a*Tl t d s aiTl r 



(3.7) 

For every equation of state p = p(n, s) the solution of this equation has the form X = 
X(IIj) = ^n = F(^/7, A 2 , {J T ) 2 )[p\, with all the dependence upon the fluid momenta in 
the quadratic form A 2 and with an extra dependence on J T = det (d r a 1 ). See Appendix B 
for the known solutions of Eq.(3.7). 

Once X = X(Ilj) is known, we have 

Y r = Y r (X,U i ) = Y r (U l ), J r = (l+n)Y r (U t )-n r .r, d T a l = -^^d r a\ 

V ~ T r (U) 

u T = u T {n,i) = —~ — , u r = u r {Ui) = 



(l + n)X(ILi) ' (l + n)X(lU) 

(3.8) 
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The energy- momentum tensor of the fluid (see Eq.(4.2) of Ref.[10]) is 



T AB (r,a)[a) = 



5S 



Vg $ 4 9ab 



(r, 3) = 



J A J B 



dn 



epU A U B -p( A g AB -eU A U B ) (r,a), 



(&„AB 



a 9cd J°J D ' 

A ttB\ 



(3.9) 



As shown in Eq.(3.11) of Ref.[2] it is convenient to express its components T TT and T T 
in terms of a mass density M. and a mass current density M. r in the following way 



ttt m = ^Hp)r) 2 -pr = ( (1+ ^ )Mi, 

T (t, o") = e(/) + P )[/ T [/ r - P V r = (- (i + n )2^ JM' ( 3 - 10 ) 

As shown in Ref.[10], by using Eqs.(3.5) and (3.6), the mass density and the mass current 
density of the fluid have the expression 



M(r,a) = ( v ^(l + n) 2 T-)(r,a) = 

> 3 g rs Y r (X,n t )Y s (X,n t ) dp(n,s) 



-(- 
-( 



X 



dn 



X \ 
+ v^p(— ,s)J(r,a) = 



g rs d r a 1 d s a 3 IL II, X X \ . 



dn 



M r (r,a) = ^(l + n) 3 g r 



T rs + n s rprr 



)(r,a)= (d r a'U t >)(T,a). 



(3.11) 



The mass current density has a universal Hamiltonian expression independent from the 
equation of state. Instead one can get the mass density Ai (r, a) explicitly in terms of the 
momenta only for those equations of state (like the dust and the photon gas) allowing to 
find an explicit solution of Eq.(3.7). Therefore in the rest of the paper we shall consider 
only the coupling of the dust to ADM tetrad gravity. 

Let us remark that Eqs.(2.8)-(2.9) imply that the energy-momentum tensor (3.9) of a 
perfect fluid looks like the one of a viscous fluid to the Eulerian observers (U A = T(w) (l A + 
m A )) 



T AB = epl A l B -pCg AB -el A l B ) + l A q B + l B q A + rr AB , (3.12) 

where p = T 2 (w) p, p = T 2 (w) p, q A = e T 2 (w) (p + p) m A and ix AB = e T 2 (w) (p+p) m A m B + 
(T 2 (w) — 1) p 4 g AB . For a viscous fluid p is the matter energy density, p the effective isotropic 
pressure (equilibrium pressure plus the associated bulk viscosity), q A the total energy flux 
vector and n AB the symmetric trace-free anisotropic stress tensor. 
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A. The Perfect Fluid Coupled to Gravity and the Bianchi Identities for the 
Energy-Momentum Tensor 

When the perfect fluid is coupled to gravity, we have to reinterpret the metric appearing 
in the action (3.14) as the 4-metric describing the gravitational field. If we add the ADM 
action to Eq.(3.14) (see the next Section for the case of dust)), then we get the Einstein 
equations 4 G AB = &2-G-T AB with the energy-momentum tensor of Eqs.(3.10) and (3.11). 
Then the Bianchi identities imply 4 V A T AB = 0. As shown in Ref.[ll], these identities may 
be written in the form 



U B [ A V A T A B ) =0, 

(5 A B -U A U B yV c T c B = 0. (3.13) 



The first line of Eqs.(3.13) implies the equation 



uA -^I + (p + P)^aU a = 0, (3.14) 



whose final form 



is obtained by using Eq.(3.5). This equation is automatically satisfied when the particles 
number conservation law A \7 A (hU A ) = and the entropy conservation law U A -J^x = are 
satisfied. 

Only three of the four equations in the second line of Eqs.(3.13) are independent: they 
can be rewritten in the form 

U B 4 V B U A = ^-U A B -U A U B )-^. (3.16) 

and turn out to be the relativistic Euler equations. They allow to express the acceleration 
a (u)A = U B A V B U A as a function of "internal forces" depending on the pressure ( for the 
dust we will get a^) A = in Eq.(4.19)) 

a ^ = Tr P {^- u ^ B )^ < 3 ' 17 > 

Finally, by using the notations of Appendix A for tetrad gravity and Eqs.(3.11) of Ref.[2], 
the identities 4 V A T AB = may be written in the following form 
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9r - E " e U »(«) A' " (1 + ») (( 3e )^ *0 + E 9 « [ E M-) M-) ^ 

+2( 3 e)- 1 E^ 3 e(a) 3 e(a)^-( 1 + n ) E 3 ^ 3e (a) 3 4) 3e (a) 3 4) T - = ' 



+ 



abrsuv 



as 

~ E ^ (H) "(«)) ( 3e ) _1 ^ - (! + «) ^ ( 3e )^ - 
as 

-E^™ T ™ + (! + ») E 9 « 3e W H) T ") 



Ef) ( 3 p, \ 3 p, s \ 3 p u 3 p v 3 p m 3 p n T 



abcuvmn 



= 0. 



(3.18) 



Since Eqs.(2.9) imply the following form of Eqs.(3.11) 



( 3 e)~ 1 M = r 2 (w) (p + p)-p, 



<*e)- 1 M r = -Y 2 {w){p + p)w r 



(3.19) 



the use of Eqs.(3.9) allows to get the following expression of Eqs.(3.18) as equations for w r 
and p 



dw r 
~d7 



dw r 



w. 



dn. 
da' 1 



: V(l + n)E^ 3 V^,-!^: + 

/ u 



+ ^(E--^ + ( 1 +»)E 



+ 



r» (p + p) 



/-, X ^P 



(dp sr-udp\ 
\8t ^ da u J 



0, 



^-E((!+-)- s +- s ) 



9 



9 (7 s 



P 



(E 3 9uvw u w v ) 



d 








T OCT 



dp 



P+(l + n)E w '^ + 



(p + p)(l + n) 3 ^-E + E' V 



w 



0, 



(3.20) 
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where 3 V is the covariant derivative inside the 3-spaces with 3-metric 3 g rs = J2 a e (a)r^( a ) 
and rf = J2 s 3 g rs n s = J2 a ri( a ) 3 e r ,, are the shift functions. 
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IV. THE DUST IN ADM TETRAD GRAVITY 



In this Section we shall consider the dust coupled to ADM tetrad gravity [1] in the globally 
hyperbolic, asymptotically Minkowskian space-times with boundary conditions killing the 
supertranslations and reducing the asymptotic symmetries to the ADM Poincare' group. 
The notations for the gravitational field are defined in Appendix A. 



A. The Action of the Dust coupled to ADM Tetrad Gravity 

Let us consider an isentro pic perfect fl uid, a dust with p = 0, s = const., and equation 
of state p = fin = ^= \J e A gAB J A J B ■ In this case the chemical potential p is the rest 
mass-energy of a fluid particle: p = mc. The full action of ADM tetrad gravity coupled to 
a dust (with positive energy) is 



S — Sadm + Sdust — p — mc, 

= Sadm- J drd 3 a m c y/e ±g AB J A J B (r, a) , (4.1) 

with Sadm given in Ref. [1]. The gravitational action depends upon 16 tetradic gravitational 
configuration variables, which are defined in Appendix A. At the Hamiltonian level there 
are 16 momenta as shown in the canonical basis of Eq.(Al). 

The dust momentum conjugate to a 1 is {T ri is defined in Eq.(3.6) with dp/dh = p) 
tt / — *\ S Sadm Y r (r,a) 

(4.2) 

The following Poisson brackets are assumed 

{a\r, a), IL(r, a')} = -5) 5 3 (a, a'). (4.3) 



As shown in Ref. [1] there are the following ten primary first class constraints involving 
only the gravitational field. 



^f(a)( T ^) ~ °> n n (r,a) « 0, 7r n(o) (r,a) w 0, 
3 M (a) (r,a) = e (a )( 6 ) (c ) 3 e {b)r (T,a) 3 n[ c) (T,a) w 0. 



(4.4) 



Moreover there are the following four secondary first class constraints 
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= 'Hgravir, 3) + M(t, (?) « 0, 



"H(a)(r, <t) = U grav{a ){T, (?) + 3 e( a) (r, <r) -M r (r, ?) w 0. 

(4.5) 



They are the super-Hamiltonian and super-momentum constraints, which depend also on 
the mass density Ai and on the mass current density Ai r = d r a 1 ITj defined in Eq.(3.11) (see 
Ref.[l] for the gravitational parts r H gra v( T ,^), Hgrav(a)( T ,&))- Eqs.(4.5) show the necessity 
of having the Hamiltonian expression of the mass density M by solving Eq.(3.7). 

For the dust we get the following solution of Eq.(3.7) (see Appendix B) 



/i(J 



T\2 



jidet 2 (d r a' 1 ) 



V vV + ^ 2 ft* det* (d r a*) + Y~ a 3 ef a) 3 e^ a) d u a* d v a> IT, IT, ' 



Y r = det(d r a i 



^ det* (d r a*) + £~~ 3 e « a) 3 e « a) du ~i dv ~i n , il, ' 



r = (l + n)y r - 3 e^n (o) J T = 



= det(d r a i ) 



(1 + n) 3 eL)^ n * 



^//i 2 cfet 2 (0 r a*) + E„^-a 3e (a) 3 ^) ^ a* S a* II, H, 



;«)]> 



J r = cfet (d r a l ). 



(4.6) 



Consequently, we can get the expression of the velocities of the Lagrangian coordinates in 
terms of the momenta by using Eq.(3.1) (whose validity in general relativity will be shown 
in Eq.(4.15)) 



d al = E, r d r a> _^ {J T E a 3 e\ a) n (B) -(1 + n) Y*) d r 



a: 



n 



(a) 



^ 2 (J T ) 2 + E~~ 3 e u (a) \ a) d u a* d v a> IT, Uj 



(4.7) 



so that the unit 4- velocity of the dust becomes ( 4 gAB U A U B = e; in the last line the 4- velocity 

is decomposed on the contra-variant ortho-normal tetrads l A (r, a), 4 E^(r, a) carried by the 
Eulerian observers as shown in Eqs.(2.8) and (2.9)) 
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u r 



J 7 



(l+n)X p(l + n)J T 

J' Ea 3e (a) 



(l+n)X p(l + n)J T 



V 2 (J T ) 2 + E 9 - ^ d " aJ ^ % 

uvija 

(l + n)X;M«) 9 - a<n *- 



- n( a ) //i 2 (J^) 2 + £ 3 ^ a) 9 U a* 9, a> Hi Uj 



uvija 



% A U A = 6 -^±^- (lp* (J T ) 2 + ^l) du a* d v <* IU n, 

V uvija 

E^Ho^4 



= -e r(iu) w r = ^ <9 r ^ n is 



^ = (i + „)tr^ + ^ 



r 4 



1 + n 



(a) 



(4.8) 



For the dust we have the following Hamiltonian expression of the energy-momentum 
tensor (3.9), of the energy density p = ep = eph and of the mass density (3.11) 



T AB = ephU A U a ^pU A U a , 



A TT B de f- TT A J T B 



M = ^(l + n) 2 T^ = (l + n) 2 ^p(Un 2 = Jl*(j T ) + E " e l) " e \a) & ^ <* j ^ n „ 
T TT M 



uvija 



P = 



M r = ^f(l + n)J2 3 9rs (T TS + n s T TT ) = d r a* U t = -p J T U r 
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T rs = pU r U s = 
e 



V 3 p r 3 p s 

l^ab e (a) e (b) 



(l+n) 2 V7 



\J^ 2 { JT ) 2 + 3 e? a) 3 e v (a) d u a* d v a> U t U 3 



(l + n)J2 3e (a) d u ol 1 n, + n (o) //i 2 (>) 2 + ^ 3 e« a) <9 U a* 9, a* II, Uj 

ui y uvija 

(l + n)J2 3 4) 5 - a ' ^ + n W /^ 2 (>) 2 + Yl 3e (a) 3e («) & & a J ' n * n : 

vj y uvija 

(4.9) 

Let us remark that for the dust, having p — 0, Eqs.(2.9) and (3.19) imply .M = 
(/>r 2 (w) p = p J T T(w) and .M r = — (/>r 2 (w) pw r = — Ai w r . 

The Hamilton equations for the dynamical matter, given in the next Subsection, imply 
the identities A V A T AB = 0, consequence of Einstein equations due to the Bianchi identities, 
whose expression was given in Eqs.(3.18) or (3.20). For the dust Eqs.(3.20) become (p = ep; 
nT = Ea 3 e( a) , 3 K rs = l 3 K 3 g rs + J2 a b a (a)(b) 3 e( a )r 3 e(b) s with the notations of Appendix 
A) 

(dw r \ - , dw r \ - dn s \ \ - ,. s dn 

TT - E " apr - E gp " (i + ») E - ^ + 

5 S / U 

+Wr (j2w u ^ + (l + n)J2™ u ™ v 3 K UV ) = 0, 



d 
d 



n; 



d 
da s 



P ~ 



(l+n)p[ ( 3 K- Y, wrwS3K ™) +J2 3V * 



w 



0. 



(4.10) 



At the Hamiltonian level they are a consequence of the Hamilton equations for the dust, 
given in Eqs.(4.14). 

B. The Dust in the York Canonical Basis 

In the York canonical basis (A2) we have (0 = 1 / 6 = (y/^y) 1 ^ 6 ) 



¥>(a) "(a) 


n n (a) 









a* 1 




7T n 7Tn( o) ~ 






n s 


nj 
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X 



4>n = (f) ^ 



yr = jr 



* ^/i 2 (J^) 2 + 0- 2 / 3 Y~ 3 Qa 2 Vua v va d u a* d v a s n, n, ' 

^ 2/3 E bsl QfVrkVskdsOtTli 



^ (J T ) 2 + 0- 2 / 3 Eauvij Qa' Vua V va O u ^ d v ^ li, II, 

r = (i + n) Y r - 3e (a) n {a) r = r <r 1/3 £ v rb 

a b 

(i + n )0-i/3Q-i^. Kb g s ^n t 

/X 2 (J T ) 2 + 0" 2 / 3 Ea««ii Qa 2 Vua Ka <9 U «9„ a* IT, 



J T = det(d r a l ), 

Tri 



1 ~ 



asuvjk 



d T a l 



E, 



^ «9 r a* _ E r ( E a M a) n w J- - (1 + n) y) d r 



a: 



rb 



/i 2 (J T ) 2 + 0~ 2 / 3 Ea,«, Qa 2 Vua Ka <9 U d v & li, II j 



M = ^ (J-) 2 + 0~ 2/3 ]T Q^Ka^a^a^a^n,, 
A^ r = ^ <9 r a 1 Hi, 



arsij 



rpr 



-5/3 



1 Qb 1 ^ra Kb 



£0 ' \ ^ 

(1 + n) 2 ^ / T^ 2 1 

"" ^ 2 (j-J +0- 2/3 Ecuvij Qc 2 V uc V vc d u a i d v am i Il j 



(i + n) r 1/3 g; 1 ^ \/ ma 9 m d n, + n (a) /^ 2 (>) + 0- 2 / 3 ^ g- 2 k c K c <9 u 9, «i n, n, 



CUVIJ 



(l + ^r^Q^Y.Vnbdn^Uj + n^Jpi^r) +0- 2 / 3 ^ g^KcKc^^^^n.n 

(4.11) 



raj 



CUVIJ 



The unit 4- velocity (4.6) and the energy density p of Eq.(4.9) of the dust are 
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U T = 



u r 



V arsij 



r 1/3 Eb Qb 1 Vrb 

H (1 + n) J T 



U T = 



- n (b) Ifj? (j t ) 2 + 0-2/3 ^ q-2 Vra Vsa dr ~i ds ~j Hl n,] , 

y arsij 

V U T + 4 g TS U s = e([(l+n) 2 -J2 4)] U T - 1/3 E Qa n (a) V sa U s ) = 

a as 
^-j^- [ L* (J T ) 2 + 0- 2 /3 Qa" Ka Ka d r «* d s IIj - 



arsij 



-*" 1 "E^«'E l '- ii -« i 4 

a si 



1 ^ Mr 



Id J' 



/i 2 (J T ) 2 + 0- 2 / 3 E arsij 2 Ka <9 r «9 S ^ IT; 11, 



(4.12) 



C. The Hamilton Equations in the York Canonical Basis 

The Hamilton equations of the gravitational field in Schwinger time gauges are generated 
by the following Dirac Hamiltonian (see Eq.(3.48) of Ref. [2] for H grav ), which depends upon 
matter only through the mass density and the mass current density 



H D = - Eadm ^ / <l rr 



/' 



T,<7 + 



+ / d 3 a 



71% - 7l( ) %(a) 

\ n n n + A" o) vr" a) ] (t,<t) = 
= H grav + J d 3 a [(1 + n) M\ (r, a) - J d 3 a ^ n r (r, a) M r (r, a) 
= H grav + J d 3 a (1 + n) m] (t, a) - 

~ j S(J E ("W ^" 1/3 Qa 1 E Ka(^)-Mr)(r )( T). 



(4.13) 
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and are explicitly given in Section IV of Ref. [2]. The expressions of the super-Hamiltonian 
and super-momentum constraints in the York canonical basis are given in Eqs. (3.45) and 
(3.42) of Ref. [2]. 

The Hamilton equations for the dust are 
d T a l (T,a) = - / dVi l + w(r,(7i) + 

+ / dVx [0-/3 £ n(a) Q-l £ Ka](r , = 
a r H > / 

= ((1 | n ) 0-2/3 ^ Qa 2 Ka Ka ^ Q* d s ^ IT 

(J r ) 2 + 0- 2 / 3 ^ 2 Kb 9 U a™ S a- n m n n 

a r 

J ' da\T,(x) 
- f dV, [#-/» 53 « W Q„- ^ V„](r, ?1 ) = 



^ Y/i 2 (^) 2 + r 2 / 3 E 6 „ Q^ 2 K fe Kb d u a™ d v a n n m n n 
/i 2 ^ T £ e '' uv e ^ k du a ' d - ak + ^ 2/3 £ 2 v ™ ^ ^ a ' n ^ Ui ~ 

uvjk asj 

~ £" 1/S £ n wG« ln *)( T '*)- ( 4 - 14 ) 



The first half of these equations allows to check the validity at the Hamiltonian level 
of the basic defining property of the Lagrangian (comoving) scalar fields a l (r, a) used as 
coordinates for the fluid flux lines 

U A (T,a) A V A a i (T,a) = U A (T,a)d A a i (T,a)=0. (4.15) 
In accord with Eq.(4.2) the inversion of the first equation gives (T ir = d r a 1 from Eq.(3.6)) 
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H J T 2/3 Erse, & Ql Vrc V sc Xj s l [d T a> - 0-V3 £ Q-l K(l d u * 

II; " — 



^(1 + nf - 0V3 Efc [g fe Esfc VshX ~l [ dr a* - ^ ^ n (c) Ke ^ 



^ (J T ) 2 + ^ q-2 Kfe y„ 6 9 U a™ d„ a™ n m U n 

^(1 + n) 2 - 2 / 3 E 6 E s , KbX^ 1 [0 T a" - n (e) Q-i Ke ^ a*] 

/iJ T (l + n) 

(4.16) 

By putting this expression for Ilj(r, a) in the second set of Hamilton equations, we 
can get the second order equations of motion for the dust variables a l (r, a): they have 
the form of three conserved currents dAjf(T,a)=0. These conservation laws imply that 
the quantities J d 3 a Ilj(r, a) are constant of the motion and generate the transformation 
a l (r,a) — > a l (r,a) + const, leaving invariant the Lagrangian (4.1). 

Finally the Hamilton-Dirac Equation (4.15) for the dust are equivalent to the 4- 
dimensional, manifestly covariant, equations 



U A *V A U B =0, (4.17) 

implying that the flux lines are geodesies. For the dust with null pressure Eq.(4.17) is 
equivalent to eq.(3.16). 



D. The 3-Orthogonal Schwinger Time Gauges 

The restriction of the Hamilton equations to the family of (non-harmonic) 3-orthogonal 
Schwinger time gauges, where 9 1 (t, a) so that the 3-metric is diagonal ( 3 g rs = <fi 2 ^ 3 Q 2 r S rs 
from Eq.(A3)), is given in Section II of Ref.[3]. The restriction to these gauges of Eqs.(4.14) 
is 
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//i 2 (j t ) 2 + r 2/3 E ^ 2 E ^ «* n *) 2 - 

y a i 



-5/3 



(1 + n) 2 I 

V^ 2 ( JT ) +0- 2/3 E c Qe 2 (E^c^n,)2 

(1 + n) r 1 / 3 Q- 1 £ <9 a a 4 n, + n (a) L 2 (>)' + ^ £ Q-2 d c a* I^) 2 ] 

i y c i 

(i + n) r 1 / 3 g,- 1 £ d fe ^ n, + n {b) L 2 (j-) 2 + 0-V3 £ q-2 ( ^ 9cQ;J - n,.)*] , 

J V c j 

tA 2 ( Jr ) 2 + ^ 2/3 E ^ 2 E ^ « l n *) 2 > 

V c i 



/i(l + 

li (1 + n) J T 



(i + n) 0- 1/3 £ Qr^r^n* 



n (r) /u 2 (^) 2 + 0~ 2 / 3 E <?a 2 E 9 " a * n ^) 2 ' 

y a i 



e { J^~ i^ 2 {JT)2 + ^ 2/3 s Q " 2 9 « al 9 « aJ n « n ^ 



a i 

f/ r = — eT(w)w r = — X] d r a l Yii, 
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d T a\T,a) = -((1 + n) 



_ 2/3 Egj Q^dga'dg^Yl, 

^//i 2 (J r ) 2 + 0- 2 / 3 £ fem „ Q 6 ~ 2 ^ « m S « n n TO n n 

a 

^ ^ Ec, ^ Qe 2 [9r ^ ~ Ea *(«) Qj 1 <9 a «*] 



IT = 

^(1 + n) 2 - 2 /3 E 6 [a E fc ^ [d T a* - Ee n {e) g-i 9e ^] 
<9 T Ui(r, a) = - E ^ ( . = 1 + U 

^//i 2 (J T ) 2 + 0- 2 / 3 £ 6m „ Q b ~ 2 ft « m ^ n m n„ 

[\ ^ JT E ^ ^ ^ ^ ^ + ^' 2/3 E 9r ^ ^ ^ - 

uvjk j 

- r^^^g-in^r^). (4.18) 

a 

E. The Hamiltonian Post-Minkowskian Linearization 

We now remember the results of Refs.[3, 4] on the Hamiltonian Post-Minkowskian lin- 
earization of ADM tetrad gravity with dynamical matter in the 3-orthogonal Schwinger time 
gauges by using the dust as matter. The asymptotic Minkowski 4-metric at spatial infinity 
is used as an asymptotic background. Post-Newtonian expansions are avoided by introducing 
a ultraviolet cutoff M c 2 for the energy of the matter. 

The basic assumption is that on each instantaneous 3-space S T we have the following 
limitation of the a-dimensional configurational tidal variables R n in the York canonical basis 

\R- a (r,a) = R (1)a (r,a)\ = 0(C) « I, 

\d u R- a (r,a)\ ~ io(C), \d u d v R- a (r,a)\ ~ -^O(C), 

\d T R- a \ = ±O(0, \d 2 T R- a \ = ^ 0(C), \d T d u R- a \ = -5,0(0, 

=S> Q a (r, a) = e S» 7- fl»(r,9) = ! + r (i) (r> &) + 0(C 2 }) 

T« = E 7a- ife, E V « ] = °> ^ = E ^ ^ ( 4 - 19 ) 

a a a 

where L is a fog enough characteristic length interpretable as the reduced wavelength \/2n 
of the resulting GW's. Therefore the tidal variables R a are slowly varying over the length L 
and times L/c. This also implies that the Riemann tensor a Rabcd, the Ricci tensor a Rab 
and the scalar 4-curvature A R behave as j^O(Q. Also the intrinsic 3-curvature scalar of the 
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instantaneous 3-spaces E T , given in Eqs.(2.5), is of order 0(C)- To simplify the notation 
we use Ra for R^)a in the rest of the paper. 

For the other gravitational variables we make the following consistent assumptions (the 
notation f\ k ^ = 0(( k ) is used) 

= 1 + 60 (1) + O(C 2 ), 
N = l + n= l + n(i) +0(C 2 ), 
n {a ) = n ( i)(a) +0{( 2 ), 

V = e(l + 2ra (1) ) +0(( 2 ), 
Vr = -en(i)(r) + C(C 2 ), 

Vs = -e 3 <7™ = -e (l + 2 (if > + 2 (1) )) 5 rs + 0(( 2 ). (4.20) 
Moreover we have (see Appendix A for the shear (7(a)(6) of the Eulerian observers) 



8ttG 



U- a (r,a) = ^n (1)B (r,*) = ±O(C)=[0 T Jk-^ 



8ti G \ -v 1 2 \ 

<7(a)(o) = ^(l)(o)(a) = 2^ 7a a n ( i )a + - 0(C ), 



<7(a)(6)|o#6 = <7(l)(a)(6)|a#> = 7 0(C), 



8ttG 



*f = 7 o(c 2 ) = £(if ) -r«) eiab a (1)(a)(fc) + -o(c 3 ), 



127rG ,„ 127rG 1 



— ^w-z 



0(C), 



4 



5 ^ 



(l)rs 



■o(C) 



= (1 - <5 rs ) <T(i)( r ) (a ) + 6 rs - 3 K {1) -d T F^ +^2(6 ra - -)<9 a 77 (1)(a) +-0(C 2 ). 



(4.21) 



For the matter we must have 
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M(r,a) = M {1) (r, a) + TZ {2) (r, a) , 

mi = M0((), J d 3 aM (1) (T,a) = McO(C), J d 3 an (2) (r, 3) = McO(( 2 ), 

M r (r,a) = M (1)r (r,a), J d 3 a M (1)r (r, a) = McO(C), (4.22) 

where M is the ultraviolet cutoff. We have M.(i)(t, <r) = ^ 0(C), where the length L is the 
wavelength of gravitational waves as shown in Ref. [3]. 

For the dust we have a l (r, a) = 0(1) and that the mass scale m = ^ must satisfy 
m = MO((). Therefore the dust momenta are first order quantities IIj(r, a) = II^^t, a) 

such that f d 3 cr M.(i) r (T,a) = f d 3 a d r a l U^j(T,a) = McO((). These conditions also 

imply Jd 3 aM {1) (T,a) = McO{Q. 

As a consequence we have the following results for the energy-momentum tensor and 
for the 4- velocity of the dust (see Refs. [3, 4] for the expression of the generators of the 
asymptotic ADM Poincare' group) 

M = M (1) + ^0(C 2 ), 

m {1) = //i 2 (^) 2 +E(E^^ n ^) 2 ' 

y a i 

M r = M {1)r = -d r a i n i + ^-0(( 2 ), 

< V; ; 11;^ H; Mc 2 

^ 2 { JT ) +Ec(Ei9c« i iii) 2 

«9 T ^ (1 )+«9 r 7W (1)r = + ^O(C 2 ), 
9 T A4 (1)r + 9 s T ( 'i s ) = 0+^O(C 2 ), 
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r 



V \ A 2 ( JT ) 2 + E a< n *) 2 + °(o, 

y c i 



U r = ^—droilli + Oit). (4.23) 

From Refs. [3, 4] we have the following solutions of the linearized equations for the gravi- 
tational field 
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0(1) ( r > °) 
n ( i)(r, <t) 



c 



c 3 A 
rAirG 1 



a 

*(l)(a)(T,*) = [| 3 ^(l) + ^^(4Ml) a -|E^Ml)c) + 

c 



2 A 

o 1 



<7(l)(o)(6)|o#& = 2 ( 9 « n (l)W + ^> w (l)(a))( r >^)- 



(4.24) 



The retarded solution for the tidal variables (the gravitational waves) and for the TT 
3-metric are (a multipolar expansion has been used; M a i = 8 a i — J2 a ^aa 7^ ^2 ji a — 

2 75fc "a); q uv ' TT is the quadrupole; see Ref.[3] for the tensors A rsuv and V rsuv ) 



4 j TT / -»\ ^ /* t3 7 3 \ * tTT / -* -* \ 

«(i)r S ( r . a ) = ~ e ~ d a l( i a 2 d rsuv {a l -a 2 ) 



T^{t-\S-& x IB 2 ) 



\a — o"i 



2G x ^ . 



n2 n uv\rr 


((r- 






(7 





+ {higher multipoles) + 0{l/r 2 ), 



R~a(r, a) = Taa rj?> (r, <?)= [r« (r, <?) = £ 7sa R- a (r, a)] 

a a 



a6 



r ( ^(T-|a- ( T 1 |,a 2 ) 



cr — OA 



o{C) 



ab 



bbuv U T H 








— * 

a 





)• 



(4.25) 



The tidal momenta can be obtained from Eq.(4.21) 
The linearization of Eqs.(4.18) is 
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d a a 1 d a a 3 Wj 



o V/ 2 ( J ) n fc) 



+ <9 a a*n ( i)( a ), 



<9 T IT = E ^ 



/i 2 J r e r ™ e ijfe <9 n a J 9, a fc + ^ d r a j IT, IU 



it = 



X ic X Jc d rOL 3 ' + 0(C 2 ), since /i = mc = McO((). 



(4.26) 



By putting the solution (4.25) for the gravitational waves into these equations we get an 
integral differential equation for the dust. 

The HPM linerization of the Bianchi identities (4.10) is (w r = Mciy/Mri) = 0(1)) 



dw r \ -v _ dw r \ -v drinks) 
~ E n ww 7w " E ^ 



dr 



1 + n (1) ) WU d - W r - 2 E ^«) 2 ^ ^(1) - E 9 r r ^ + 

u « s 

+ ( E - u ^ + a + ») E ^ 3 %)™ ) - ^ 



d ^ d 

n (!)( s ) 



da s 



p = 0. 



(4.27) 



Let us remark that in the non relativistic limit U A >->■ (l,^( r nr )), w r >->■ —U^ nr ^ and 
P | — >■ P(nr), where p( nr ) is the non relativistic mass density. Therefore, in the non relativistic 
limit, the previous equations become the non relativistic Euler's equations (with n(i) being 
the Newtonian gravitational potential) and the mass conservation equation 



du U £^ 0tff„) 



dn 



(i) 



dr +2^ U (nr) da . 



(nr) da s da r ' 

P(nr) = 0. 



(4.28) 
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V. KINEMATIC AL AND DYNAMICAL ASPECTS OF THE DUST 

In this Section we study the congruence of the dust flux lines, its connection with the 
skew congruence associated with the 3+1 splitting and the Eulerian point of view for the 
dust. Then we will define the acceleration, the expansion, the shear and the vorticity of the 
dust flux lines. Finally we will face the problem of how to select the subset of the irrotational 
motions of the dust and we will study the connection between the resulting dust 3-spaces 
and the 3-spaces of the Eulerian observers of the 3+1 splitting. 



A. The Congruence of Dust Flux Lines 

As said in Section IIC there are two congruences of time-like observers associated with 
each 3+1 splitting of space-time, i.e. with every global non-inertial frame. Their unit 

4-velocities are P(t,<t) = (*£Z A )(T,a) and v»(r,a) = -^L=(t,o) = (z%v A> )(T,a). 

When the dust (or every type of perfect fluid) is present we also have the congruence of 
the time-like flux lines with unit time-like 4- velocity U^(t, a) = z^(r,a) U A (r,a), which in 
general is not surface-forming. 

The flux line of the dust emanating from a Q at r = are denoted 

cju>(a) = Z>1 ( T W = ^(W A ) ; = c ( W A ))- (s.i) 

Here A is an affine parameter, C(off )(^) = ( r (^)5 ^( r (^)> <?o) ) > r (0) = 0, S(0, & ) = a , 

C(o a t o )(0) = (O'^o)- This equation defines the 3- coordinates S r (r, a ) identifying the location 
at time r of the flux line emanating from a r at r — 0. 

The flux lines are the integral curves of the unit time-like 4-velocity, namely they are the 
solutions of the equations 



= U A (r = r(A) = C ( V- o) (A); ° r = Z r (r(\),a ) = C(W A ))' 



(5.2) 



If we would impose the three gauge fixings on the gravitational field 

U r (r,a) « 0, 

^ (4.12), 

n (a) (r,a) « (r,a)), 

y/l* (J T ) 2 + 0- 2 / 3 Earsij Qa" Ka Ka d r ^ d s ofl IT, Uj 

a\T,&) « a* (a), from £g.(4.15), (5.3) 
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we would get U^{t, a) = t> M (r, a), namely the fluid 4-velocity would coincide with the 4- 
velocity of the observers of the skew congruence of the 3+1 splitting, whose world-line 
passing through a r Q is now also a flux line with constant 3-coordinate a r Q 

4M X )) = C5, ffo) (A) = ^(r(A),a ), 

=► C[ 0A) (A) = 57(r(A),^) = (7S. (5.4) 

Since the conditions (5.3) determine the inertial shift functions, it means that there is a 
choice of 3-coordinates 9 l (r,a) on the 3-spaces S r implying this coincidence 8 . 

For r = A Eqs. (5.3) are consistent with the first half of the Hamilton equations (4.14). 
In this case we have S r (r, a ) = a r and a l (r, a) = P{&). 

The comoving coordinates of the gauges (5.3) are usually used in cosmology when the 
matter is dust. See for instance Ref.[17] for the dust and Refs. [7, 18] for more general 
fluids. 

B. Other Aspects of the Dust Flux Lines 

As shown in Ref. [11], instead of the (comoving) Lagrangian coordinates a l (r, a) describing 
the fluid in the 3-spaces S T , we can use the set of Eulerian coordinates S r (r, a ) appearing 
in the definition (5.1) of the fluid flux lines. 

If we identify the afnne parameter A with the time r of the observer origin of the radar 
4-coordinates, i.e. if we put A = r, the flux lines are described by the radar 4-coordinates 

CU)(r)=(r; £(r,<7 )). 

Since Eq.(4.15), i.e. U A (r,a) <9a cc 1 (t, <t) = 0, imply that the scalar fields a l (r, a) can be 
interpreted as labels assigned to the flux lines (they are constant along the flux lines), we 
can write 

a l (r,S(r )( T )) =a\0,a o ) = f al(a ). (5.5) 

The quantities a. % (a ) are Cauchy data on the Cauchy surface S r=0 at r = 0, whose 3- 
coordinates are denoted a r Q = S r (0,(? o ). If we invert the relation a r = S r (r, a a ) to get 
°~o = 9%(t,<t), we have 

c/(r, a) = a% g s (r, a)) = a*(^g(r, <?)). (5.6) 
Since the conserved particle number is Af = f v ^ ^ d 3 o~ J T (a. l {r, a)), see after Eq.(3.2), 
on S r=0 the initial number density is 

n (a ) = r(a%ff)) = det (^^)- (5.7) 
If the fluid has compact support we can take a*(a ) = a r =l , so that n (a ) = 1. 



As said in Appendix A, the r-preservation of the three gauge fixing for the angles 9 l implies three equations 
for the determination of the shift functions. 
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C. The Eulerian Canonical Coordinates 

The flux lines of Eq.(5.1), Q^- ^(r) = ^r, E(r, <x )j with A = r, are identified by the 
3-coordinates o" r = E r (r, <x ), where are the 3-coordinates on the Cauchy surface E r= o. 

As shown in Ref.[ll], there is a point canonical transformation allowing to pass from the 
(comoving) Lagrangian canonical coordinates a l (r, a), IIj(r, a), to a set of Eulerian canonical 
coordinates S r (r, a D ), K r (T,a ). Therefore the new fluid coordinates coincide with the 3- 
coordinates a r Q (and not with a {a )) on the Cauchy surface S r=0 and for r > they label 
the flux lines with E r (r, c? ) instead that with a l (r, a) = a (g^(T,(x )). The initial number 
density h (a ) is assumed to be a known function given as part of the Cauchy data. Also 
ot l (a ) is known from the Cauchy data. 

The relation between the old coordinates and the new ones is obtained by rewriting the 
relation a l (r, <?) = a % (g^(r, a )) in the following form 

« l (r, a) = J d 3 a Q det (r, a ) 5 3 (a r - TT{r, <j )) a* (<f ). (5.8) 

By using a*(r, E(r, a )) = a (<J ), it can be shown that the inverse of Eq.(5.8) is 

E r (r, a ) = J So det (r, a) 5\a {a ) - a\r, a)) a r . (5.9) 

The generating functional and the new momenta of the point canonical transformation 

are 



37,11,] = Jd'aU^a) Jd 3 a det[^j (r,a ) 5 3 (a r - E r (r,a )) a (a ), 



K r (r,a ) = 



5$ 



= " det (^) ^ff^^f^ /^ n ^-) 5 V-S s (r,a )). 

(5.10) 

Instead for the inverse canonical transformation we have the following generating func- 
tional 
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a\ K r 



d 3 a K r (T,a) / d 3 <x det 



r, <r) <5 3 (a*(a ) - a l (r, a))a r , 



U t (r,a) 



4 



5$' 



<5a; 1 (r, a) 

_det (la) ^'^la 7 ^'^ J d3(T oK r ( y r,a )5 3 ( y ai( y a ) -a(r,a)). 



(5.11) 



It can be checked that these expressions for the momenta are one the inverse of the other. 

By using the notations introduced after Eq.(3.6), i.e. Z ir (r,<j) = <9 r a*(r, <r), we get 
V I^(t,S(t,<7 )) 8sr fl (r f o) = ^v 21 with «o(^o) defined in Eq.(5.5). Since we have = 



rfoUgo) . . da* (T,Y,(T,a )) _ 8tt'(T,S(r ; ff„)) 

(<9 T a* = — ^ r yfT ir ) the following results 



£ r X, r (r,£(r,a ))^M we get from Eq.(3.2) 



J-(r,S(r,a )) 



9a'(r , S(r , a )) J r (r,E(r,a )) 



dr 



J-(T,E(T,a )Y 



det (d r a%r,J:(T,a )) = -det" 1 — - (r,a ) det ' " l 



J r (r,S(r,a )) 



-Q-^j( T i°o), => J T (r, a) = [h(a ) det 

„ , /9S«\ 057(r,* o ) 

n (a )det | — ((r.ff,,) 



9 as 
! /<9£ r 



9 as 



vo=gs(T,a) 



dr 



(5.12) 



As shown in Ref. [11], these results imply the following form of the mass density of the 
dust 



M(t, a) 



det" 



0*2 



/i 2 n 2 (a ) + V'fo S(r, <r )) K r (r, a ) K s (r, a Q ) 



J d 3 a 5 3 (a r - Z r (r,a )) ^ h 2 (a ) + V s (^, S(r, <?„)) K r (r, a ) K s (r, a ). 



(5.13) 



Instead for the mass current density we get the following result valid for any kind of 
perfect fluid (for the dust we have M. r = —fj, J T U r ) 
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U r {T,a) = 5.14 
lin {g^(T,a)) 

D. Acceleration, Expansion, Shear and Vorticity of the Dust 

The covariant derivative of the dust 4- velocity U (t, a), whose Hamiltonian expression 
is given in Eqs.(4.12), allows to find the associated acceleration a^, expansion 0(u), shear 
C(u)ab and vorticity u)(u)ab of the dust. The analogous quantities for the Eulerian observers, 
associated with the 3+1 splitting and having the unit normal l A (r,a) to the 3-spaces E T , 
are given in Eqs.(A5)-(AlO) of Appendix A. 

For the dust we get 

4 VacUb = eUAa(u)B + V(u)ab + ^9(u) 3 h(u)AB — W(u)ab, 3 'h(u)AB = 4 9ab — e UaUb , 
a fu) = U B 4 Vs U A =0, a,(u)A = A 9ab afu), afu) Ua = 0, 

&{u)ab = V(u)ba = — 2 i a (u)A Ub + a,{u)B Ua) + - ( 4 Va Ub + 4 Vb Ua) — - 0(u) 3 h(u)AB=, 
= 2 ( 4 Va Ub + 4 Vb Ua) — - 0(u) 3 h(u)AB, <?(u)ab U b = 0, 

W(U)AB — —U(U)BA = VABCD ^(u) U° = 

= ~\ (a<u) A Ub - a (u)B U A ) - | ( 4 V A U B - 4 V B U A )=-\ ( 4 V A U B - 4 V B U A ), 
u (h) = 2 f] ABCD U(u)bc Ud, U(u)ab U b = 0, Ua = 0. 

(5.15) 

We have explicitly shown the consequences of the fact that the equations of motion (4.17) 
imply the vanishing of the acceleration a^(r, a)=0. 

In many contexts, especially in cosmology, one considers irrotational dust, having 
uJ(u)ab( t ,&) — 0. In this case the unit time-like 4-velocity U A (r,a) is surface forming, 
namely there is a 3+1 splitting of space-time whose 3-spaces £(t/) T are orthogonal to U A (r, a). 
In these cases there should exist gauges of canonical gravity in which the congruence of Eu- 
lerian observers coincides with the dust congruence of observers. 
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However our dynamical dust has in general non- vanishing vorticity. If we try to impose the 
condition oo(u)ab{ t i <?) = by asking to have Ua(t, B) ~ Ia(t, B) — e ^1 + n(r, B)j (1; 0), we 
get the conditions U r (r, B) 0. But from Eqs.(4.13) this implies £V <9 r a l (r, a) IIj(r, <r) ps 0, 
i.e. Ilj(r, a) ps since J T (r, a) 7^ 0. Then the second of Eqs.(4.14) gives J2 r &r ({} + 

n) e ruv tijk d u a 3 d v a k j (t, a) ~ 0, i.e. three conditions on the lapse function implying n = 
n(r). 

To understand better this restriction, in the next Subsection we will explore how to 
implement the condition u)(u)ab(t,B) = as a restriction on the space of solution of the 
dust motion. Then, after this identification, we will explore which gauge-fixing on canonical 
gravity will imply U A l A for the irrotational motions of the dust. 

E. The Subset of Irrotational Dust Motions 

From Eqs.(5.15) we get that the conditions W(u)ab( t , B) ~ are implied by W(i/) rs (r, B) 
0, because we have U(u) Tr {r, B) = J2 S (^{u)rs jf^j ( r ) B) ~ 0. 

Therefore we must study uj(u)rs= — | ( 4 V r U s — 4 V S C/ r ) = — f (9 r U s — d s U r ) d = — | fi(i7) rs ~ 
as a restriction on the solutions for the dust motion. 

Let us remark that this restriction is preserved in time by the Hamilton equa- 
tions, because Eqs.(4.13), (4.14) and (4.12) imply ^^^={Q {u)rs {T,B), H D } = 
[-gf? (Q(u)su jp) + gf? (tyi/)™ (t,^) ~ 0. Therefore it is enough to impose the con- 
dition fi(£/) rs (0, B) ^ on the Cauchy surface S r= o- This also shows that these conditions 
can be interpreted as first-class constraints to be added by hand to select the family of 
irrotational motions of the fluid. 

Let us remark that £l(jj) rs ~ are only two independent conditions, because we have the 
identity d r £l{u)uv + d u Q(u)vr + d v £l{u) ru — 0. 

If we make the following non-covariant decomposition of U r (r,B) 9 (A = ^ r the 
function c(r) is required by the boundary conditions at spatial infinity and will make the 
gauge fixing later introduced in Eq.(5.25) well defined) 



Compare with the analogous decomposition of the electro-magnetic vector potential in Refs. [10], [3], 
where it is shown that one must firstly fix the 3-coordinates on the 3-space S T and then make the non- 
covariant decomposition. 
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U r (r, a) 



= d r S(r,a) + U± r (r,a), 



S(r, a) = c(r) + £ | C/ r (r, a) = c(r) - i £ | ^^(r, a), 



d 

U± r (r,a) = ^n {u)vr (r,a), ^ d r U ±r (T,a) = 0, (5.16) 

-y r 

the two conditions of vanishing vorticity are 

U ±r (r,a)ttO, U r (r,a)^d r S(r,a). (5.17) 

These are the two first-class constraints to be added by hand to eliminate the states 
of motion of the dust with non-zero vorticity. They show that two of the three pairs of 
canonical variables a l (r, <?), IL;(t, a), describing the dust are associated with vorticity. The 
addition of two suitable gauge fixing constraints (whose form is not known), so to get two 
pairs of second class constraints, would allow to go to Dirac brackets and to describe the 
irrotational dust with only one pair of canonical variables. 

Since S(t, a) and U± r (r, a) have non-trivial Poisson brackets (i.e. {S(t, <?), S(r, a')} ^ 0, 
{^(u)rs(r, a),Q(u)uv(i~, <?')} 7^ 0) {S(t,v),Q{u)uv(i~, <?')} 7^ 0), it is not easy to identify the 
two gauge fixings to be added. This is not easy also in the Eulerian point of view, where 

Eqs.(5.14) and (5.12) imply U r (r,a) = ^ T ^^'^ . The presence of the initial number 

density and of the function gg(r, a) imply a complicated expression for U± r (r,a), so that 
the two gauge fixings will be complicated functions of S r (r, a ). 

However it is possible to identify the reduced phase space of the irrotational dust also 
without knowing the gauge fixings explicitly. Eqs.(4.12) imply 

m M{ r, 3 -)) = - y^jj) (r. W - O, 
{det(S„a'(T,?)),S(T,9')} = --i 3 (<r" -"'"). 

{det(d r a\r, «?)), n (u)uv (r, a')} = 0. (5.18) 
As a consequence we get 

{U r (r,a),U s (r,a')}\ n(u)rs=0 = 0, 
{S(r,a),S(r,a')}\ nmra=0 = 0, 

{S(r,a),n (u)uv (r,ff')}\ n(u)rs=0 = 0. (5.19) 
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Even if we do not know the expression of the two gauge fixings, it turns out that the 
reduced phase space of the irrotational dust is spanned by the two conjugate canonical 
variables 

u(r,a) =h J T (r,a) = 11 det(9 r a i (r,a)), S(r, a) = ^ d r U r (r, a), (5.20) 

r 

because, due to Eqs.(5.18) and (5.19), their Dirac brackets are 



{z/(r, a), i/(r, a')}* = {u(t, a), i/(r, v')}n {u) „=o = 0, 

{S(r, a), S(r, a')}* = {S(r, a), S(r, a')}n {u)rs =o = 0, 

Mr, a), S(r, a')}* = {u(r, a), S(r, a')} Q(u)r3=0 = -5 3 (a r -a' r ). 



(5.21) 



While z/(r, a) describes the numerical density of the dust (M — - J d 3 a z/(r, <?)), the func- 
tion S(t, a) allows to put the dust unit 4- velocity in the form ( 3 g rs = 0~ 2 / 3 J2 a Qa 2 Vra V sa , 

n " = Ea Ha) Qa 1 V ra ) 



U T (r,a) = (1 + n(r, a)) ^1 + 3 g rs (r, a) d r S(r, a) d s S(r, a) + n r (r, a) d r S(r, a), 



U r (T,a) = d r S(r,a), U a q Ub = 1. 



(5.22) 



With the substitutions J T = det(<9 r a: J ) >->■ ^ u, d r a l Hi \-> —vd r S, the part f d 3 a ^(1 + 
) M + n r M r ^j (r, a) of the Dirac Hamiltonian (4.13) takes the form 



H. 



(R) 



dust 



d 3 a v(t, a) 



1 + n(r, a)) ^1 + 3 g rs (r, a) d r S(r, a) d s S(r, a) + n r (r, B) d r S(r, a) 



(5.23) 



This reduced Hamiltonian gives the following Hamilton equations in the reduced phase 
space of the irrotational dust 

8 T S(t, ff)={S(r, a), H^J* = ((1 + n) ^/l + 3 g™ d r S 8 S S + n r d r s) (r, a) = U T (r, a), 



v (1 + n) 



3 g rs d s S 



y/l + 3 g rs d r Sd s S 



n 



(5.24) 



which allow to write Ua(t, g)=8a S(t, a) consistently with Eqs.(5.22). The second half of 
Eqs.(5.24) allow to check the constancy of the particle number: d T Af — - j d 3 ad T v(t, er)=0. 
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The result Ua(t, <?)=<9a S(r, a) shows that the unit time-like 4- velocity of the irrotational 
dust is orthogonal to the 3-spaces S(t, a) = const. Moreover one can check the validity of 
Eq.(4.17) in the reduced phase space, namely that the flux lines are geodesies. 

If we add the constraint 

X (r,a) = S(r,a)-S(r)^0, i.e. c(r) « S(r), (5.25) 

implying U r (r,a) ps 0, Ua{t,o) p=* Ia{t,o) and IIj(r, a) pa for the irrotational dust in the 
original phase space, its preservation in time produces the constraint 

(R) -| * dS(r) _ TT t _ _, N dS(r) 



d T X (r, a) = { X (r, a), H^Y ~ = U T (r, a) 



dr dr 

pa l + n(r,a) « 0, 

dr 

^ (5.24) 

<9 T u(t, o) pa — ^ d r {vrf^j (r, <r), n = n(r). (5.26) 

r 

But this means that there is a gauge fixing for the inertial gauge variable 3 K(t, a) of canon- 
ical gravity whose r-preservation implies that the lapse is only a function of time. As a 
consequence there is a 3+1 splitting of space-time whose 3-spaces coincide with the ones of 
the irrotational dust. 

If before the restriction to irrotational dust we had chosen the 3-coordinates 6 l (r, a) such 
that Eqs.(5.3) hold, i.e. U r (r, a) ~ 0, then Eqs.(5.26) imply fi( a )(r, a) ~ and d T z/(r, er) ~ 0. 

This, with the extra condition n(r) = 0, is often the starting point in cosmology with 
irrotational dust and comoving 3-coordinates [7, 17]. In this formulation the irrotational 
dust is not described by canonical coordinates like a l (r, a), IIj(r, a), but only by a function 
p (r, a) such that the energy-momentum tensor is T^ v = p £/ M U u . This function satisfies 
the Bianchi identities T^ v '.„ = 0, which imply 10 d T [p 0] = (0 = y/ del 3 g) and the result 
d T p = 3 K p . 

In our approach we have the function p of Eqs.(4.12), whose Hamilton equations are 
induced by Eqs.(4.14) and by the Hamilton equations for the gravitational field [2]. By 
using the Hamiltonian Bianchi identities of Ref.[2] and by restricting to irrotational dust 
satisfying both U r ps and U r pa (i.e. with n pa pa 0) we recover d T p = 3 K p. 

Finally let us remark that by using Eqs.(5.15) we could also study: i) dust without 
expansion: it requires the study of the equation 4 V U a ~ 0; ii) shear-free dust: one should 
study the equations f ( 4 Va U b + 4 Vb Ua) ~ § 0(u) 3 ^(i/)ab- 



10 This has the consequence that the mass inside an averaging volume is constant My(r) = 
J v d?<j[p$\{T,a)=M v {T ). 
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F. About Irrotational Motions in Arbitrary Perfect Fluids 



Let us make some remarks on the existence of irrotational motions for arbitrary perfect 
fluids like the one discussed in Appendix C. 

In Ref. [19] it is shown that Eqs.(5.15) together with the definition ^ 4 V j 4 4 Vb — 

4 Vb 4 V j 4^ Uc = —^RcdabU d of the 4-Riemann tensor imply the following equations for 
acceleration, expansion, shear and vorticity i\AB] means anti-symmetrization) 

3 h {u)A C3 h {u)B D (u E V E U(u)CD - 4v [ c a ([/)D ]) + 
e 2 

+2(T (U)e[a W V) B ] + 2 9 V)U(U)AB = 0, 

(5.27) 



Let us look whether a perfect fluid with arbitrary equation of state may admit irrotational 
motions, namely motions with null vorticity 

W(u)ab = 0. (5.28) 

While the second of Eqs.(5.27) is identically satisfied, the first of Eqs.(5.27) implies that 
the neccessary condition for the existence of irrotational motions is 

3 h { u)A C3 h { u)B C4 V^a {u)D ^ =0. (5.29) 

Therefore Eq.(5.29) is a restriction on the equation of state of the fluid, giving the pressure as 
a function of the energy density. Moreover Eq.(5.29) must be compatible with the relativistic 
Euler equations (see Eqs. (3.16) and (3.17)). 

This problem was considered in Ref. [18], where it was shown that if for a family of motions 
the 4- velocity of the fluid admits the following parametrization in terms of a scalar function 
S(r, a) 

U A (r,a) = —!—V A S(T,a), 
h{T, a) 

h(r, a) = vV B (r,<?) V A S(r,a) V B S(r,a), (5.30) 
then the acceleration turns out to depend only on the normalization function h(r, a) 



a(u)B(r,a) = U A (r, <?) V a U b (t, a) = 

= (V - U b (t, B) U A (r, a)) V A In h(r, a), (5.31) 
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and the condition (5.29) is satisfied. For the dust Eqs.(5.22) and (5.24) show that Eq.(5.30) 
is satisfied with h(r,a) = 1. 

As shown in Ref. [18] barotropic fluids, with equation of state p = p(p), admit a family of 
irrotational motions for which Eqs.(5.30) and (5.31), and then (5.29), are valid (the function 
h being proportional to U(p) = j -fe). 
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VI. CONCLUSION 



Brown's formulation [9] of perfect fluids allows to describe them only in terms of three 
Lagrangian (comoving) coordinates in Minkowski space-time by means of an action principle 
whose Lagrangian is determined by the equation of state of the fluid. This action was 
reformulated as a parametrized Minkowski theory in Ref.[10] and then studied in the rest- 
frame instant form of dynamics. The main drawback of this approach is that we can get an 
explicit closed form of the Hamiltonian quantities only for few physically relevant equations 
of state, including the dust and the photon gas, due to the necessity of solving a trascendental 
equation (see Appendix B). 

As a consequence, in this paper we studied the coupling of the dust to ADM tetrad gravity 
in globally hyperbolic, asymptotically Minkowskian space-times. We found the Hamiltonian 
formulation of the gravitational field with dynamical (not test) dust as matter and we gave 
the Hamilton equations of dust in the York canonical basis of Refs.[l, 2]. In this way 
we can disentangle the inertial gauge effects of the gravitational field from the tidal ones 
(the gravitational waves in the linearized theory). Also the Hamiltonian Post-Minkowskian 
linearization [3, 4], avoiding the Post-Newtonian expansion, is given in the 3-orthogonal 
Schwinger time gauges. 

By using radar coordinates adapted to a time-like observer we define the Hamiltonian 
theory in global non-inertial frames, centered on the observer, with well-defined instanta- 
neous 3-spaces, dynamically determined by Einstein's equations. In the York canonical basis 
the basic inertial canonical gauge variables are: i) three angles describing the freedom in 
the choice of the 3-coordinates inside the 3-spaces (their gauge fixing determines the shift 
functions); ii) the York time (the trace of the extrinsic curvature of the 3-spaces) describing 
the general relativistic remnant of the freedom in clock synchronization (its gauge fixing 
determines the lapse function). 

In this framework we have studied the problem of selecting the subset of the irrotational 
motions of the dust. For this subset there are special 3-spaces determined by the dust 4- 
velocity and we have studied the problem of which gauge fixing is needed for having these 
dust 3-spaces coinciding with the 3-spaces of the global non-inertial frame. Also the Eulerian 
point of view [11] was discussed in the case of dust. 

Since dust is the type of matter used in cosmological models (usually irrotational dust in 
comoving coordinates) and since there are indications that at least part of dark matter can 
be explained as a relativistic inertial effect induced by the York time [4, 20], the material 
of this paper is preparatory for the ADM formulation of cosmology. In particular we want 
to use this Hamiltonian description of dust in the framework of back-reaction [7, 21], in 
which an averaging procedure inside the 3-spaces (leading to a breaking of homogeneity 
and isotropy) opens the possibility to describe dark energy, and all the associated effects 
as the accelerated expansion of the universe, as an effect induced by the non-linearities of 
Einstein's equations. Since the spatial average is well defined only for 3-scalar functions 
and since most of the quantities describing the gravitational field and the dust are 3-scalars 
in the York canonical basis, we now have a framework for studying the spatial average of 
most of the Hamilton equations. Moreover we can study which cosmological notions have 
an inertial origin being functionals of the inertial gauge variable York time. These problems 
will be faced in a future paper. 

Moreover we have to study the photon gas in the framework of this paper. Also we can 



44 



look for Hamiltonian expressions of the fluid mass density, whose Lagrangian expression 
gives approximations to the equations of state usually used for fluids (an inverse procedure 
to the one used in Appendix B). This would open the possibility of studying compact fluid 
bodies without symmetries (relevant for the description of star) and their multipoles [22, 23] 
at the Hamiltonian level and to see whether our approach can be useful in formulating a 
well-posed Cauchy problem with a free boundary for the body, a still unsolved problem at 
the mathematical level [24]. 
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Appendix A: Canonical ADM Tetrad Gravity and of the York Canonical Basis 

In this Appendix we review the formulation of canonical ADM tetrad gravity developed 
In Refs. [1] starting from the ADM action considered as a functional of the cotetrads defined 
in Eqs.(2.1). Then we introduce the York canonical basis and the expansion and the shear 
of the congruence of the Eulerian observers. 

1. The Original Canonical Variables of ADM Tetrad Gravity 

As said in ref. [1-3] and with the notations of Eqs. (2.1)-(2.5), in ADM canonical tetrad 
gravity the 16 configuration variables are: the 3 boost variables <P( a )(r, a); the lapse and 
shift functions n(r,a) and 7i( a )(r, a); the cotriads 3 e( a ) r (r, a). Their conjugate momenta are 
^ip (a) (T, a), 7T n (r, a), 7r n(o) (r, <?), 3 tt^(t, a). There are 14 first-class constraints: A) the 10 pri- 
mary constraints of Eqs. (4. 4): 7i> (a) (Y, a) 0, 7r n (r, a) « 0, 7r n (r, a) ~ and the 3 rotation 
constraints M( a )(r, a) ps implying the gauge nature of the 3 Euler angles a( )(r, a); B) the 4 
secondary super-Hamiltonian and super- momentum constraints / H(r,a) ~ 0, W(a)(j, a) m 
of Eqs. (4.5). As a consequence there are 14 gauge variables (the inertial effects) and two 
pairs of canonically conjugate physical degrees of freedom (the tidal effects, which become 
the gravitational waves in the linearized theory). 

The basis of canonical variables for this formulation of tetrad gravity, naturally adapted to 
7 of the 14 first-class constraints, is (see Refs, [1, 25] for the direction- independent boundary 
conditions at spatial infinity) 



<P(a) 


n 


n (a ) 


3 e {a )rl 




« 


Ti n W 


n (a) 


3 -u 1 



{n(r, a), ir n (r, a)} = 5 3 (a, a), 

{n (o )(r, a),n n(b) (T,a)} = 5( a )( b) 5 3 {a,a), 

{<P(a)(T,a),n nb) (T,a)} = 5 {am 5 3 (a,a), 

{ 3 e (a)r (r,a), 3 7r ( s b) (r,<T')} = 5 {am 5 s r 5 3 (a, a). (Al) 
2. The York Canonical Basis 

In Ref. [1] we studied a point canonical transformation on the canonical variables (Al), 
implementing the York map and adapted to 10 primary first-class constraints (it is a Shan- 
mugadhasan canonical transformation adapted also to the constraints M( )(r, a) ~ 0). It is 
based on the fact that the 3-metric 3 g rs is a real symmetric 3x3 matrix, which may be diag- 
onalized with an orthogonal matrix V(6 r ), V' 1 = V T (J2 U V ua V ub = 5 ah , J2 a V uaV va = 5 UV , 
^2 UV e wuv V ua V vb = J2 C e abcV cw ), detV = 1, depending on 3 Euler angles 6 r n . The gauge 



Due to the positive signature of the 3-metric, we define the matrix V with the following indices: V ru . 
Since the choice of Shanmugadhasan canonical bases breaks manifest covariance, we will use the notation 
V ua = J2v Ki« &v(a) instead of V u ( a y We use the following types of indices: a = 1, 2, 3 and 5 = 1,2. 
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Euler angles 6 r give a description of the 3-coordinate systems on S T from a local point of 
view, because they give the orientation of the tangents to the 3 coordinate lines through each 
point (their conjugate momenta are determined by the super-momentum constraints). We 
only consider 3-metrics with 3 distinct positive eigenvalues. In this way we get the following 
York canonical basis 



<P(a) 


n 


n {a) 


3 e( a )r| 




s 


7T n P 


a 


n n (a) ' 


« 


3 *u 1 



¥>(a) "(a) 


ra n (a) 






^1 




VT n ~ 7T S(a) « 






n„| 



(A2) 



In the York canonical basis we have (from now on we will use V ra for V ra (9 n ) 

def 

to simplify the notation; we use the following definitions: ri( ) = ^ b ^(a)(6)( a (c)) 

3 e( a )r d = Et %)(6)(«(c)) 3 e(6)r, 3 e( a) = J Efc R (a)(b) («( c )) 3 e[ 6) , where i? (a)(b )(a; (c )) are rotation 
matrices, i? T = 



4 

9rr 

4 

9rs 




c[(l + n) 2 -X;n(«) 

a 

-£ ^(a) 3 e( a )r = ~C 1/3 QaV ra n {a) , 

a a 

-e 3 g rs = -e 2 / 3 ]T Q 2 V ra V sa , Q a = ^ *», 

a 

6 = ^ = Jd^g~ = % 3 e (a)r = 1 / 3 Q a V ra , % } = r 1 / 3 Q- 1 V ra , 



n U = %)w( Q; (e)) 7r (b)> 

6 

0" 1/3 [v; a g- 1 (0^ + E^ n s) + 



7T, 



(a) 



+ YY QT 1 Vrl€ahVwt ; B tw Trf } 



7T, 



7T7 = 



Cmir &(a)l "(a) ' 

imra 

„3 1 

Yhx~G ~ 3^e ^ 6(a)r ' 



Ha = 7aa 3 vr[ a) 3 e (a)r . 



(A3) 



47 



The set of numerical parameters ^- aa satisfies [1, 13] J2 U 7s« = 0, J2 U 7a«75« = 
Yla laulav = $uv — §■ Each solution of these equations defines a different York canoni- 
cal basis. See Ref. [1] for the SO(3) Cartan matrices Aij(9 n ) and B = A -1 . 

In Eq.(A3) the quantity 3 K(t, a) is the trace of the extrinsic curvature 3 K rs (r, a) of the 
instantaneous 3-spaces £ r , whose expression in the York canonical basis is 

a f, 
ab twi ^Wa Va<V b 

The York internal extrinsic time 3 K(r, a) is the only gauge variable among the momenta: 
this is a reflex of the Lorentz signature of space-time, because ir^ and 9 n can be used as a set 
of 4-coordinates [1]. Its conjugate variable, to be determined by the super-hamiltonian con- 
straint, is = 6 = 3 e, which is proportional to Misner's internal intrinsic time; moreover 
4> is the volume density on S r : Vr = J R d 3 a(p 6 , R C S T . Since we have 3 g rs = 4> 2 ^ 33 g rs with 
det 3 g rs = 1, is also called the conformal factor of the 3- metric. The two pairs of 3-scalar 
canonical variables n s , a — 1,2, describe the generalized tidal effects, namely the inde- 
pendent degrees of freedom of the gravitational field. In particular the configuration tidal 
variables R & depend only on the eigenvalues of the 3-metric. They are Dirac observables 
only with respect to the Hamiltonian gauge transformations generated by 10 of the 14 first 
class constraints. 

~ (0) 

Since the variables <p and tx\ are determined by the super-Hamiltonian and super- 
momentum constraints, the arbitrary gauge variables are a( a ), (p( a ), & l , tt^, n and n^ a y As 
shown in Refs.[l], they describe the following generalized inertial effects: 

a) a.{a)(j, a) and </?( a ) (r, it) are the 6 configuration variables parametrizing the 0(3,1) 
gauge freedom in the choice of the tetrads in the tangent plane to each point of S r and 
describe the arbitrariness in the choice of a tetrad to be associated to a time-like observer, 
whose world-line goes through the point (r,a). They fix the unit ^-velocity of the observer 
and the conventions for the orientation of gyroscopes and their transport along the world-line 
of the observer. 

b) 9 l (r, a) describe the arbitrariness in the choice of the 3-coordinates in the instantaneous 
3-spaces S r of the chosen non-inertial frame centered on an arbitrary time-like observer. 
Their choice will induce a pattern of relativistic inertial forces for the gravitational field , 
whose potentials are the functions V ra {9 1 ) present in the weak ADM energy Eadm [1, 2]. 

c) fi( )(r, a), the shift functions appearing in the Dirac Hamiltonian, describe which points 
on different instantaneous 3-spaces have the same numerical value of the 3-coordinates. They 
are the inertial potentials describing the effects of the non- vanishing off-diagonal components 
A g rr {r,a) of the 4-metric, namely they are the gravito-magnetic potentials responsible of 
effects like the dragging of inertial frames (Lens-Thirring effect) in the post-Newtonian 
approximation. The shift functions are determined by the r-preservation of the gauge fixings 
determining the gauge variables 6* l (r, a). 

d) 7t^(t, a), i.e. the York time 3 K(t, a), describes the non-dynamical arbitrariness in 
the choice of the convention for the synchronization of distant clocks which remains in the 
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transition from special to general relativity. The York time (being a momentum) gives rise 
to a semi-definite negative kinetic term in the Dirac Hamiltonian [1-3]. This unusual inertial 
effect is connected to the problem of the relativistic non-dynamical freedom in the choice of 
the instantaneous 3-space, which has no non-relativistic analogue (in Galilei space-time time 
is absolute and there is an absolute notion of Euclidean 3-space). Its effects are completely 
unexplored. 

e) n(r, a), the lapse function appearing in the Dirac Hamiltonian, describes the arbitrari- 
ness in the choice of the unit of proper time in each point of the simultaneity surfaces E r , 
namely how these surfaces are packed in the 3+1 splitting. The lapse function is determined 
by the r-preservation of the gauge fixing for the gauge variable 3 K(t, a). 

3. The Expansion and the Shear of the Eulerian Observers. 

Let us now consider the geometrical interpretation of the extrinsic curvature 3 K rs of the 
instantaneous 3-spaces E r in terms of the properties of the surface-forming (i.e. irrotational) 
congruence of Eulerian (non geodesic) time-like observers, whose world-lines have the tangent 
unit 4- velocity equal to the unit normal orthogonal to the instantaneous 3-spaces E T . If we 
use radar 4-coordinates, the covariant unit normal el a — (1 + n) (1; 0) of Eqs.(2.4) has the 
following covariant derivative 

4 Va e Ib = cIa 3 ob + o ab + - h AB — u AB — el a 3 a B + 3 Kab, 



3 K AB = 3 K rs b r A b%, b A — S r A + 3 e( a ) ri( a ) S A , h A B = 4 9ab — e I a Ib- 

(A5) 

The quantities appearing in Eqs.(A5) are: 
a) the acceleration of the Eulerian observers 

a —l Vjji — g OB-, &A — «r^, 

3 „ 3 -r 3 -u „ 3 „ 3 -r 



a r = d r In (1 + n) = a r e\ a) e u {a) n u = a r e\ a) n (a) , 

'(a) e (a) a s, 



3 o r = -e 3 e[ a) 3 ef a) 3 a s , V = 0; (A6) 



b) the vorticity or twist (a measure of the rotation of the nearby world-lines infinitesimally 
surrounding the given one), which is vanishing because the congruence is surface- forming 

10 AB = -ujba = 7T {I a 3 a B - Ib 3 oa) - 77 ( 4 Va Ib ~ 4 V B I a) = 0, 



ojab l B = 0, u A = ^- 4 rj ABCD ubc h = 0; (A7) 
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c) the expansion 12 , which coincides with the York external time, is proportional in cos- 
mology to the Hubble parameter H 13 and determines the dimensionless (cosmological) de- 
celeration parameter q = 3 l A 4 Va | — 1 = — 3 9~ 2 l A 8a — 1, 

4 A 3 4,GV^_ 12nG 
= VaI =-e K = - 3 r _ -e — — vr^, 

H = \ e =\ lA ^Al = -e^ n , q = 3l A *V A 1 --l; (A8) 



d) the shear 
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^ab = o~ BA = -^ ( 3 a A Z B + 3 a B I A ) + I ( 4 V A Ib + 4 V B Za) - ^6 3 h AB = 

= ( 3 K rs - 1 3 <? rs 3 iT) SI S B , V* <t ab = 0, o~ab l B = 0. (A9) 

By explicit calculation we get the following components of the shear along the tetrads 
(2.2) 



oja) o(/9) ( o_A _B 

&AB = CT(a)(/9) #A -^B = fl'AC 9B»^ , °"(a)(/3) = °~ AB #( a ) E^, 

°"(o)(o) = 0, a (o)(a) = 0, 

C(a)(6) = 0"(6)(a) = (^rs - ~ Vs 3 -K") 3 C( a) 3 e( b) , ^ (7( a )( a ) = 0. (A10) 

a 

<7( a )(6) depends upon 6 r , 0, i? a , Trr 6 *-* and II a . 

As a consequence, by using Eqs.(A3) and (A4) we have 



It measures the average expansion of the infinitesimally nearby world-lines surrounding a given world-line 
in the congruence. 

I is a representative length along the integral curves of 4 E^ , describing the volume expansion (contraction) 
behavior of the congruence. 

It measures how an initial sphere in the tangent space to the given world-line, which is Lie-transported 
along the world-line tangent (i.e. it has zero Lie derivative with respect to Z M d^), is distorted towards 
an ellipsoid with principal axes given by the eigenvectors of g^ v , with rate given by the eigenvalues of 
<r"„. 
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7 8irG ^ ^ c3 iv^ 

F(«)(a) = — Taalla, -> = ^ ^aa V(a)(a) , 

a a 

7 I 87T G sr-^ €abt V wt \r-^ ($) 

*" wwU = " t- -o. or 1 V ' ' 

c 3 

=>• 7r i 6,) = t^jt^ ^ ^ ^ a ^& 1 e * afc °"( a )( 6 ) l a ^ 6 ' 

3^ 3 3 3 / ~\ 

K rs = — - g rs # + o"( a )(6) e( a ) r e(6) s . (All) 

Therefore the diagonal elements of the shear of the Eulerian observers describe the tidal 
momenta n s , while the non-diagonal elements determine the variables Ttf\ determined by 
the super-momentum constraints. Moreover their expansion 9 is the inertial gauge variable 
determining the non-dynamical part (general relativistic gauge freedom in clock synchro- 
nization) of the shape of the instantaneous 3-spaces E T . 
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Appendix B: The Fluid Velocity in terms of the Fluid Momentum 



Let us restrict the action (3.4) to isentropic, s = const., perfect fluids. As shown in 
Ref. [10] the Lagrangian and the fluid momenta have the form 15 



n = 



\J\ def X \J\ def f( X 

— p = p( n )=p( ^ ,_ ) = ///(—), 



(l + n)^7 yfy' 



(l + n)^7 <fy 



X 



L = -n(l + n)yfyf(—), 
U l (r,a) = ^-(r,a) = -p 



df(x) 



dx ,x= v^ X 



(Bl) 



Eq.(3.7) becomes 



X 2 



df(x) x , 
dx x= ~7f' 



= // 2 ( J T )\ A 2 = Y, V s d r a 1 n, d s a' Uj, 



(3.11) 



X = ^n = F(^,A\(rY)[p}, 
A 2 X 



M = 



2 df(x) l 
dx \x= 



rsij 



(B2) 



Some possible equations of state for barotropic fluids, i.e. with p = p(p), are (in the 
isentropic case one gets p = p{n) by solving p(p(n)) = n^^- — p(h); the definition of the 
sound velocity is v 2 = c 2 ^y-)' 



1) p — 0, dust: this implies 



p{n) = pn = p , i.e. 



/(-) = — 
'V a/7' 



dX 



1 

7' 



(B3) 



The equation for X and its solution are 



15 In general ji is not the chemical potential but only a parameter. The quantities X, Y r and Tu are defined 
in bq.^.bj. We nave w |at = i + „ - (i+ n )x - (i+n)x > 9^ ~ ~ J (i+n)x ■ 
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X 2 [p 2 (J T ) 2 + A 2 ] = B 2 , 



x = 



^ 2 (J T ) 2 + Ersij V s d r a* d s ai IT, IT, ' 
\J T \ Zsi V'dsa'IU 



^2 ( jry + 3 g u V du al dv ~j Ut n . ' 

M = L 2 (J T ) 2 + ^ 3 g rs d r a' d s a? IT Ii r (B4) 



2) p = kp{n) = n^^- —pin) [k ^ — 1 because otherwise p = const., p = 0). The previous 
differential equation for p{n) implies 

P(n) = ph k+1 = p{^-) k+ \ i.e. 
V7 



[for fe-^Owe recover case 1)]. More in general one can have k = k(s): this is a non-isentropic 
perfect fluid with p = p{n, s). 

The equation for X is 

^-Vl^ + i^T^ 1 ^ 2 (B6) 

In general this equation cannot be solved explicitly, but is soluble for k = 
|, |, 1, 2, — |, — |, — 2: the equation for X is linear for k = 1, quadratic for /c = 2, |, cu- 
bic for k — |, — |, — 2 and biquadratic for A; = — |. 

By using Section V of Ref.[10] we get the following results for the cases k = 1, 2, |, | (A 2 
is defined in Eq.(A2)): 

2a) k = 1, p = p, p = pn 2 



x = — (^) 2 + ^ 



2pJ r y 7 

2 7 -V 2 (J-) 2 -( 7 + i)A 2 
2/iv/7(^) 2 



= " ^.'t J 2 ' . (B7) 
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2b) k = 2, p = 2p, p = ph 3 



X 



V2\ 



i + Wi- 



4 7 2 A 2 



7 V 



X 2 + 



9p 2 (J T ) 6 ' 
7 2 A 2 



3/i 2 (J-) 2 1/ 



(B8) 



2c) fc = \, p = \ p, p = pn 3/2 \ 



X 



Kf 2(jT)2 + 



7 A 4 7 1 / 2 A 2 



/i vy^ v 7 9/i 2 (j-) 4 3/i(j-) 2 7' 

2 y/2 ^ 



2d) photon gas, k — h, P = h p, p = ph^ 3 : 



x-\r\\ 3A2 ^ 1/3 1 (a 27A ^ i 27 A% 7 V /3 

1 '[ 16 (J T ) 8 / 3 2 1 / 3 V 2 1 V 6 (^ T ) 8 V 2 1 V 6 (^ T ) 8 / 
~ 2173 (~ 



( 27 ^ 6 7 / 27 A 6 7 y/3 



2 1 V 6 (^ T ) 8 V 2 1 V(-A) 8 / 

^ = n- 1/6 i 2/3 (i 2 ^|g). (Bio) 



3) p = kp" 1 [n) = n^rp^- — p(n) (7 7^ 1) [26]. It is an isentropic polytropic perfect fluid 
(7 = 1 + i). The differential equation for p(n) implies [a is an integration constant; the 



chemical potential is // = || 

ah an 
pin) = — = j — , i.e. 

[1 - kian)"!- 1 }— [l-fc(ara)»]» 

/(4=) ~ 



y/1 [1 - k(aj=p-^ [l-k(a^~' 



¥ = > - *<^r* - ^[1 - *(-^)*]-^. (bid 

Instead in Ref.[27, 28] a polytropic perfect fluid is defined by the equation of state 
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k(s) 

p(h, n s) — p h H (p n) 7 , 

7 — 1 



(B12) 



and has pressure p = fe(s)(/in) 7 = (7 — l)(p — fin). The (in general non explicitly soluble) 
equation for X is 



X 2 



/ Y 1 \ 2(n+l) 



B 2 



(B13) 



In Ref. [10] there is also a discussion of the relativistic ideal non-isentropic (Boltzmann) 
gas [29] (p = nk B T, p = mc 2 nT{(3) - p, p = *±£ = mc 2 T (/?)). 
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